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!> ■ ABSTRACT 

^ I We investigate the possibility of substantial inflation of short-period Jupiter-mass 

' planets, as a result of their internal tidal dissipation associated with the synchronization 

. and circularization of their orbits. We employ the simplest prescription based on 

■ an equilibrium model with a constant lag angle for all components of the tide. We 
, show that 1) in the low-eccentricity limit, the synchronization of the planets' spin 

I with their mean motion is established before tidal dissipation can significantly modify 

"fH I their internal structure. 2) But, above a critical eccentricity, which is a function of 

Qh| the planets' semimajor axis, tidal dissipation of energy during the circularization 

O ' process can induce planets to inflate in size before their eccentricity is damped. 3) For 

• moderate eccentricities, the planets adjust to stable thermal equilibria in which the rate 

■ of their internal tidal dissipation is balanced by the enhanced radiative flux associated 
J> . with their enlarged radii. 4) For sufficiently large eccentricities, the planets swell 

^ , beyond two Jupiter radii and their internal degeneracy is partially lifted. Thereafter, 

^ I their thermal equilibria become unstable and they undergo runaway inflation until 

" " ' their radii exceed the Roche radius. 5) We determine the necessary and sufficient 

condition for this tidal inflation instability. 6) These results are applied to study 
short-period planets. We show that for young Jupiter-mass planets, with a period less 
than 3 days, an initial radius about 2 Jupiter radii, and an orbital eccentricity greater 
than 0.2, the energy dissipated during the circularization of their orbits is sufficiently 
intense and protracted to inflate their sizes up to their Roche radii. 7) We estimate 
the mass loss rate, the asymptotic planetary masses, and the semi-major axes for 
various planetary initial orbital parameters. The possibility of gas overflow through 
both inner (LI) and outer (L2) Lagrangian points for the planets with short periods 
or large eccentricities is discussed. 8) Planets with more modest eccentricity (< 0.3) 
and semi-major axis (> 0.03 — 0.04 AU) lose mass via Roche-lobe overflow mostly 
through the inner Lagrangian (LI) point. Due to the conservation of total angular 
momentum, these mass-losing planets migrate outwards, such that their semi-major 



- 2 - 



axes and Roche radii increase while their mass, eccentricity, and tidal dissipation rate 

decrease until the mass loss is quenched. 9) Based on these results, we suggest that the 
combined effects of self-regulated mass loss and tidally driven orbital evolution may be 
responsible for the apparent lack of giant planets with ultra-short periods ^ 3 days. 
10) Mass loss during their orbital circularization may also have caused the planets 
with periods in the range ^ 3 — 7 days to be less massive than long-period planets 
which are not affected by the star-planet tidal interaction. 11) The accretion of the 
short-period planets' tidal debris can also lead to the surface-layer contamination and 
metallicity enhancement of their host stars. 12) Among the planets with periods of 
1-3 weeks today, some may have migrated outwards and attained circular orbits while 
others may have preserved their initial eccentricity and semimajor axis. Therefore, 
planets with circular orbits are expected to coexist with those with eccentric orbits 
in this period range. 13) Gross tidal inflation of planets occurs on the time scale 
~ 10® yrs after their formation for a brief interval of 10^ yrs. The relatively large 
sizes of their classical and weak-line T Tauri host stars increases the planets' transit 
probability. The inflated sizes of the tidally heated planets also increases the eclipse 
depth of such transit events. Thus, the tidal inflation and disruption of planets may 
be directly observable around classical and weak-line T Tauri stars. 

1. Introduction 

Recent searches for planets around nearby solar type stars, based on the radial velocity 
method, show a roughly logarithmic period distribution in the currently observed period range 
of 3 days to several years (Marcy et al. 2000). The absence of ultra-short-period Jupiter-mass 
planets with P < 3 days, which corresponds to an orbital semi-major axis a of 0.04 AU, appears 
to be real and not an observational selection effect (Gumming et al. 1999). The host stars of any 
ultra-short-period planets are expected to have large radial velocity modulations which would 
make them conspicuous. Glose-in planets also have greater transit probabilities. The unsuccessful 
photometric search for transiting planets around main sequence stars in the core of the globular 
cluster 47 Tuc (Gilliland et al. 2001) provides more evidence for a lack of Jovian planets with 
ultra-short orbital periods. Thus, the cutoff in the extrasolar planets' period distribution is likely 
to be caused by some physical effect. For example, it has been suggested (Kuchner & Lecar 2002) 
that the inner disk might be hot enough to allow the magneto-rotational instability to operate and 
thereby is largely depleted, terminating the inward migration through planet-disk interaction and 
causing a lack of giant planets with orbital periods less than 3 days. In the case of 47 Tuc, most 
of the Jovian planets might have been expelled from their host extrasolar systems due to frequent 
stellar encounters in the dense stellar environment (Davies & Sigurdsson 2001), or they may have 
had difficulties in forming in circumstellar disks with low metallicities (Pollack et al. 1996). 

In this contribution, we examine the effect of tidal interaction between ultra-short-period 
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planets and their host stars and explore the possibihty that this physical process may have led 
to the apparent cutoff in the observed period distribution. This investigation is an extension of 
our earlier calculations of the interior structures of weakly eccentric Jovian planets at constant 
orbital distances under the influence of interior heating and stellar irradiation (Bodenhcimcr, Lin, 
&; Mardling 2001; hereafter BLM). In these previous calculations, the interior heating is imposed 
to be constant in time and is uniformly distributed within the planet. BLM suggest that this 
heating flux may be extracted from the tidal dissipation associated with the synchronization of 
the planets' spin with their orbital mean motion and the circularization of eccentric orbits with 
small semi-major axes a ~ 0.04 AU. Under these assumptions, BLM show that Jovian planets can 
be inflated to equilibrium sizes larger than two Jupiter radii, given an adequate rate of interior 
heating. In the proximity of their host stars, the Roche radii of the short-period planets do not 
greatly exceed their physical sizes. Thus, gross inflation of any planet could lead to tidal disruption 
and Roche lobe overflow. We investigate this process and discuss its implications. 

Mass loss induced by Roche lobe overflow may also be important for the formation of 
short-period planets. Shortly after the discovery of the short-period planet around 51 Peg (Mayor 

& Queloz 1995), we suggested that it may have been formed at a distance comparable to the 
orbital radius of Jupiter (several AU's), after which it migrated toward its host star (Lin ct al. 
1996) as a consequence of disk-planet interaction (Goldreich &: Tremaine 1980; Lin & Papaloizou 
1986; Takeuchi et al. 1996). The survival of the short-period planets requires the termination of 
their orbital evolution. Among many scenarios. Trilling et al. (1998) suggest that the size of 10^ 
yr old Jupiter-mass protoplancts is more than twice that of Jupiter, so that they would overflow 
their Roche lobe if their orbits decayed to within ~ 0.03 AU from their host stars. Trilling et al. 
(1998) further postulate that the lost material flows through the LI point and angular momentum 
is transferred from it to the planet's orbit, which would balance the angular momentum loss to 
the disk and halt the orbital decay. 

In §2, we briefly recapitulate the time scale and energy dissipation rate associated with both 
synchronization and circularization of planetary orbits. We describe the response of the planetary 
envelope in §3 and describe the mass-loss flow pattern associated with Roche-lobe overflow in 
§4. Based on some analytic approximations, we illustrate the dynamical evolution of tidally 
unstable planets in terms of a sequence of stages in §5. We also present in §6 the results of 
numerical calculations of various models of synchronization and circularization processes. Finally, 
we summarize the results and discuss their implications in §7. 



-4- 



2. Tidal interaction between short-period planets and their host stars 

The equation of relative motion between a short-period planet with a mass Mp and its host 
star with a mass M* and separation r is 



+ (1) 



1=1 

where Mt = Mp + M^,, and the acceleration fj includes the relativistic potential of the host star, 
the contribution from other planets, and the tidal and spin distortions of the star and the planet 
(Mardling & Lin 2002). 



2.1. Synchronization of the planet's spin 

For computational simplicity, we focus our attention in this paper on the evolution of the 
spin angular frequency (f^p,* for the planet and star, respectively), e, and a of one single planet 
in the limit that its spin and orbital axes are both parallel with the spin axis of the star. We also 
assume the rotational and tidal distortion of the star and the planet is approximately symmetric 
about the line separating them so that they do not significantly affect the orbital evolution. The 
only remaining contribution which can lead to angular momentum transfer and induce orbital 
evolution is the acceleration due to the tidal damping within the planet and the star where 

where n = {GMt/a^y/'^ is the mean motion of the mutual orbit of the planet and the star, and 
-Rp,*) ^p.*5 and Qp^tf are respectively, the size, Love number, and the Q-value, of the star (subscript 
*), and the planet (subscript p) (Eggleton et al. 1998). The above expression is based on the 
derivation for equilibrium tides with a constant tidal lag angle. 

The planet's specific orbital angular momentum h = r x r changes at a rate 

^...f. (3) 

where ft = fp + f>^. The corresponding total rate of angular momentum transfer from the spin of 
the planet/star to the planet's orbit 

where /j, = M^Mp/M^ is the reduced mass, and Ip^^ = ap^^ep^^Mp^^^Rp ^ is the moment of inertia 
of the planet and the star, respectively. The quantity ep^^, represents the fraction of the planet 
and star which participates in the tidally induced angular momentum exchange. Since gaseous 
planets have extensive convection zones which can be fully mixed, ~ L But, solar type stars 
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have shallow surface convection zones so that e* ^ a few times 10~^. The coefficients 0;^,* are 
determined by the internal structures of the planet and the star. The equation of state appropriate 
for Jupiter-mass planets can be approximated with a polytrope with an index 1 (Hubbard 1984), 
in which case, Up = 0.26. 

If the planet's/star's spin is aligned with the orbit, we find from equations (2) and (4) that 

^v,*={ . n. l (7^] f%^V[/3(e)n-/4(e)Op,.]+u;,,. (5) 



P'* y (l-e2)i/2M, jy a J { n J .^^^^^ { MpMt ) \ 2g;,,(l - e^)!/^ 
where 



(6) 



/3(e) = (l + fe^ + fe^ + ^e«)(l-eV, (7) 

/4(e) = (l + 3e^ + ^e^)(l-eV/^ (8) 

Q'p^, = 'iQp,*/'^kp^^: is the modified Q-value (cf. Murray &: Dermott 1999), 
Jo = /x|h| = Mp{GMta{l — e^))^/^ is the total angular momentum of the planet's orbit, 
iOp{= —Ip^p/Ip) is the changing rate of planetary spin due to its structural adjustment, and oj* 
is the changing rate of spin due to angular momentum transfer, by processes such as the stellar 
winds, between the stars and the external medium (Dobbs-Dixon, Lin, &; Mardling 2002, hereafter 
DLM). Although the contribution from ujp is generally negligible, that from cj* due to the stellar 
wind can be important for the evolution of fi* and a. In the limit of small e, the rate of change of 
the planet's angular frequency is 

which induces the planet to evolve to a state of synchronization with Qp = n. 



2.2. Circularization of the planet's orbit 

Tidal friction also induces the Runge-Lenz vector, e = r x h/GMt — r, to change (Mardling 
& Lin 2002) at a rate 

de _ /2(ffr)r-(r-r)ft-(ffr)r 



dt V GM, J- ^''^ 

The magnitude and direction of e correspond to the planet's orbital eccentricity and longitude of 
periapse so that when averaged over an orbit, the rate of eccentricity change reduces to 

de 

-^=9p + 9* (11) 
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where 

/i(e) = (l + 15eV4 + 15eV8 + 5eV64)/(l -e2)i3/2^ and /2(e) = {I + ^e^ /2 + /^) / {I - e'^ f . In 
the limit of small e and nearly synchronous planetary spin, the above expression can be further 
reduced (Goldreich &; Soter 1966) to 



de 
dt 




(13) 



In the above equation, we only considered the planet's contribution in f. When the dissipation 
of the planets' tide in the star is also taken into account, the damping rate is modified from the 
expression in equation (13) by a factor 

For planets with relatively large Rp, /? ~ 1 which indicates that the star does not significantly 
modify the e damping rate. If the host star has the same values of Q' and average internal density 
as those of the planet, its contribution would be weaker than that of the planet by a factor of 
Rp/R^. However, the contribution from a very rapidly spinning host star, with 

can reduce the total rate of e damping or even induce a net rate of e excitation with a positive 
de/dt (DLM). 

Equations (2), (4), and (10) indicate that the evolutionary time scales for Q and e are 
proportional to the magnitude of Q'p ^ . These expressions are derived using equilibrium tidal 
models with constant lag angles. But, in reality, the stellar tidal perturbation is likely to excite 
gravity and inertial waves. These waves which propagate within the planet are dissipated 
either through convective turbulence in its interior or radiative damping near the surface. The 
dissipation of the waves also leads to the deposition of angular momentum which in general leads 
to differential rotation. The nature of these complex processes are poorly understood and the 
theoretical estimates of the Qp,* of Jupiter-mass planets and solar-type host stars remain uncertain 
(Zahn 1977, 1989; Goldreich & Nicholson 1977, 1989; loannou k Lindzen 1993a,b, 1994; Lubow 
et al. 1997; Terquem et al. 1998; Goodman &; Oh 1997). Some aspects of these uncertainties 
are associated with how and which waves are excited and dissipated. Using a prescription for 
dissipation due to convective eddies, Zahn (1989) found an expression for the tidal torque which 
when scaled with the expression in equation (2) gives a frequency dependence for 

Q' =^Mp!^(^Y (16) 
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where tf = [ME?/L]^^^ is the tidal dissipation time for a convective body, and A2 ~ 210~^ when 
it is fully convective. But, using a slightly different prescription for the dissipation in convective 
eddies, Terquem et al. (1998) find an expression for the tidal torque which corresponds to a Q'p 
which is independent of n. Even less certain is how the interior of the gaseous planets and stars 
may readjust after the angular momentum deposition may have introduced a differential rotation 
in their interior. 

The theoretically determined (5*-values using both approaches imply a slow rate for close 
binary stars to circularize their orbits, which is inconsistent with the observations of solar-type 
binary stars in stellar clusters of various ages (Mathieu 1994). Mathieu et al (1992) also pointed 
out that a power-law function (e/eP^^/^) provides a close fit to the slope of the observed cutoff 
period. Such an empirical expression would correspond to a Q' value which is independent of n. 
Nevertheless, the observations are equally well fitted with an index of 10/3, and an index of 16/3 
cannot be ruled out. 

To carry out a comprehensive analysis on the physical process associated with the dissipation 
of the tidal disturbance is beyond the scope of the present paper. In view of these uncertainties, 
we took the simplest prescription based on an equilibrium tidal model with a constant lag angle 
for all components of the tide. This may not necessarily be the most reliable model. All the 
uncertainties associated with the physical processes are contained in the Q' values and we shall 
parametrize our results in terms of it. For a fiducial value, we note that the Q^-value inferred for 
Jupiter from lo's orbital evolution is 5 x 10"^ < Qp < 2 x 10^ (Yoder & Peale 1981). With this 
Qp-value, orbits of planets with Mp and Rp comparable to those of Jupiter, and with a period 
less than a week, are circularized within the main-sequence life span of solar-type stars. If we use 
equation (16) to extrapolate from the observationally determined values for Jupiter, Q!p ~ 10^ 
which would imply that only systems with periods less than 3 days would be circularized within a 
few Gyr. But, if all planets are formed in situ or brought to their present location shortly after 
they are formed, their extended radii would enhance the rate of tidal circularization. Within the 
observational uncertainties, all known extrasolar planets with periods less than a week and some 
with periods less than two weeks have circular orbits. 

With the theoretical estimate of stellar Q^-values (~ 10®), the time scale for orbital evolution 
of Jupiter-mass planets is longer than the stellar age of a few Gyr unless their orbital period is less 
than ^ 8 — 9 hours (Rasio et al. 1996, Marcy et al. 2000). But, using the observed circularization 

time scales as a calibration, the inferred value of Q'^ is ~ 1.5 X 10^ for the very young binary stars 
(Lin et al. 1996) and 10^ for the mature binary stars (Terquem et al. 1998). In this paper, we 
shall adopt the observationally inferred values for „, and discuss the implications of our results 
if they are somewhat larger. 
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2.3. Tidal evolutionary time scales 



With these values of Qp,*, the synchronization, circularization, and migration of the planet's 
orbit occurs on time scales (tq, Te, Ta) which are much longer than the orbital period, and they 
can be obtained from the time-averaged equation of motion. In the limit of /3 ~ 1, small e and 
nearly synchronous planetary spin. 




(MA 

\Mj) 



lday/(27r/np) 



O.34(M*/M0)V2(O.O4AU/a)3/2 _ lday/(27r/np) 



m7 



0.04 Auy 



Rj_ 

Rn 



Myr, 



0.33 



(Mp\ (Mq 



3/2 



13/2 



(17) 



(18) 



\de/dt\ ' ) \Mj J \MJ V0.04AU 

The above derivation for Xg is based on the approximations that the planet's spin is nearly 
synchronized with its mean motion and the tide raised on the planet by the star provides the 
dominant contribution to the e evolution. The magnitude of Tg would be shortened/lengthened 
by the planet's tide raised on a slowly /rapidly (relative to Opc) spinning star (DLM). Prom these 
equations we find 



lap€p 




Op 


2 




fip — n 



(19) 



In the limits that e << 1 and Mp « M*, r^p << Tg, and the nearly-synchronous approximation 
in equation (18) is justified self consistently. 

Under the action of tidal torque and the conservation of total angular momentum, Jo = Jp + J^. 
In the case without mass transfer between the star and the planet, this exchange also leads to an 
evolution of the semi-major axis at a rate 



a = a 



2ee , 2(J^ + Jp 



(1 - e2) 



+ 



Jo 



(20) 



where is the usual planetary migration time scale through disk-planet interaction. Since the 
planet quickly establishes a state of synchronous rotation, the magnitude of Jp becomes much 
smaller than that of J* and Tg = |e/e| << tq^ < \ Jo/J*\ < \Jo/Jp\ so that a declines together with 
e. But, after the orbit is circularized, a and J* + Jp evolve together with the same sign and the 
evolution of a is mainly driven by the dominant effect of the tide raised on the star by the planet 
(Goldreich & Sotcr 1966, Mardling & Lin 2002) such that 



_ a QiM^ {a/R 





a — . — 




13/2 



f 1 Gyr 

-ft* 



(21) 

in the limit = oo. Provided Q'^ ~ Q'p and /3 ~ 1, Tq > Te > Tf^. Por i?* = Rq, Ta is much larger 
than Td (except for a < 0.01 AU), and is also much larger than the relevant time scales tr and Tg 
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for planetary inflation and eccentricity damping (which are ~ 1 Myr, see the discussions in §3.3, 

§5, and the simulation results in §6.2) that concern us in this paper. Thus, we can neglect the 
evolution of a due to the tide raised on the star by the planet today. But, during the classical T 
Tauri phase when i?* ~ 3Rq , is a significant fraction of the breakup angular frequency, and 
~ 1.5 X 10^, the magnitude of Tq is comparable to that of rd(~ 10^~^ yr) at 

such that protoplanets' migration induced by their tidal interaction with their protoplanetary 
disks may be stalled by that with their host stars (Lin et al. 1996). Although the magnitude of 
Td may differ from its assumed value of ^ 10^ yr in the above estimate, ahait has a very weak 
dependence on r^. The value of ahait is also a weak function of Q'^. Two orders of magnitude 
increase in Q'^ would reduce Ohait by a factor of < 3. 



a = Ohalt - 



3. Planet's Structural Adjustment due to its Internal Tidal Energy Dissipation 

In this section, we consider 1) the effect of tidal dissipation in increasing the internal 
temperature of the planet, 2) the necessary conditions for planetary inflation, 3) the sufficient 
condition for Roche lobe overflow. 



3.1. Energy budget 

During the synchronization and circularization processes, kinetic energy is continually 
dissipated into heat. From equations (20) and (6) we find that both e and lead to the evolution 
in a and a corresponding change in the orbital energy Eorb = —GM^Mp/2a at a rate 

-E'orb = — -E'orb-- (23) 

a 

In addition, the planet's spin energy -Egpin = evolves at a rate 

4pin = 4pin + E%:^ = apMpRpRpnl + IpUptlp (24) 

where the first term on the right hand side represents the planet's size adjustment (see below), 
and the second term corresponds to the changes in the kinetic energy. Prom equations (6) and 
(20), we find 

Etide = Eorh + E^pin = f^n^a^ (t^) + ^P^P^P = ^(^ - ^^P X r)-ft (25) 

is the total tidal energy dissipation rate within the planet in its rest frame (Eggleton et al. 1998, 
Mardling &; Lin 2002). (A similar expression can be derived for the rate of energy dissipation 
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within the star, but its effect on tfie stellar structure is negligible.) Prom equations (5), (11), (12), 
and (25), we find 

^''"^ = - (^) (^) (^) ' H^'^'^ - 2nn,h,{e) + n'h,{e)] , (26) 

where /13(e) = (1 + Se^ + 3e^/8)(l - e^)-^/^ /14(e) = (1 + 15eV2 + 45e^/8 + 5eVl6)(l - e^)-^, 
and /15(e) = (1 + 31eV2 + 255e^/8 + IBSe^/lG + 25e^/64)(l - e2)-i5/2. 



For small e's, the leading terms of e^ and Qp/n — 1 in equation (26) are reduced to 



v 

2^3 



i^tide - I ^ 



— - 1 + e + 23 

n J \ Zn'^ n j 



(27) 



The first term in the square bracket on the right hand side of the above equation represents the 
energy dissipation due to the synchronization of the planet's spin. If planets are formed as rapid 
rotators, similar to Jupiter and Saturn, equation (27) reduces to 

E^^,:^E,^^J-^^^h^^. (28) 

When synchronization is established the terms inside the last bracket are reduced to le^ /2 and 
equation (27) may be approximated as (c.f. Goldreich & Soter 1966) 

£;tide-£^te = 29 

are 

We now consider the effect of tidal dissipation on the total energy budget of the system. In 
response to this source of energy, the planet adjusts its internal structure to evolve toward an 
energy equilibrium in which the rate of change of the planet's total energy is given by 

-E'tot = -E'bind + -^spin = — (30) 

where L is the rate of loss of energy from the system via radiation from the planet's surface. In 
this derivation, we assume that the planet maintains a state of hydrostatic equilibrium. The tidal 
energy dissipation rate is contained in E'tot in the following manner. The binding energy has three 
contributing components 

GMl 

^'^^'^'^ = -"^""rM ~ '^'Mt) ^ ^^^^ 

where q^^p = 3/(5 — n^^p) and n*^p is the polytropic index for the star and planet respectively. For 
a polytrope of index 1, q^^p = 3/4 and for a uniform distribution of density, g*_p = 3/5. 

Neglecting any changes in the gravitational energy of the host star, 

_qpGMlR^ GMpM^a 

Rp{t) Rp^ 2a{t) a ^'^ ' 
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Note that Mp in the above equation can also be a function of time when the size of planet is 
inflated up to its Roche radius and therefore mass loss occurs. Similarly, the spin energy can be 
separated into two pieces (see eq. 24). In principle, rotation may also affect the binding energy of 
the planet. But from equations (31) and (24), we find the energy ratio 

where 

is the planet's Roche radius and Mj and Rj are Jupiter's mass and radius, respectively. After 
the planet attains synchronization within its Roche radius, its rotational energy is relatively small 
compared with its binding energy. Since -Egpin << -E-bind; the effects of Coriolis and centrifugal 
force are negligible in the calculation of the planet's internal structure. In the above expression for 

we introduce a factor of (1 — e) to represent the tidal radius at periastron where the stellar 
tidal effect is at an maximum. 



3.2. Tidal inflation instability 

Tidal dissipation within the planet leads to changes in a and Op, and we have already defined 
the rate -Etide in equation (25). This process also leads to the internal heating of the planet at a 
rate — E'tide- Prom the above equations, we find that the planet's size would change at a rate 

^ = qpGM^/Rl + apMpRpnl ^^^^ 

provided the adjustment would not directly modify Op (Mardling &; Lin 2002). 
In a thermal equilibrium, Rp = and equation (35) reduces to 

C = Cg = — £^tide (36) 

such that the surface radiative flux is balanced by internal energy dissipation. In response to a 
small perturbation 6Rp from its equilibrium radius 

Re = Rp{£ = — -Etide) (37) 

and in the case of Tg >> tr (i.e. e can be roughly treated as a constant with time for this small 
perturbation), a planet adjusts its radius at a rate 



R„ 



/9//na2\ ( MA /Re\^ ( n^e) - ^nnph^e) + n^h^ieU (SR. 



\ 2Q'p ) \Mp) \a) \ \pGMllRl^apMpRp^l rJ ^^^^ 
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where 

Note that Rp and 5Rp would have opposite signs and the thermal equilibrium would be stable if 
7 > 5. But if 7 < 5, small increases in Rp would lead to runaway inflation. 

In general C \s a, function of Mp, Rp, and the existence of the core. Based on three sets of 
equilibrium models, BLM obtained an approximate Re — C relations in which 

log|^ = A{Mp) + B{Mp)\og^ + C{Mp) Uog-^) ' (40) 

it0 L,Q V LqJ 

where, in the range lO-^L© < £ < lO-^L©, A = (3.11, 1.012, -0.269), B = (1.011,0.4689,0.1631), 
and C = (0.0642,0.0289,0.00978) for Mp = (0.63,l,8)Mj respectively. In the high luminosity 

limit, i?e becomes a more rapidly increasing function of C (with a decreasing 7 magnitude) 
as degeneracy and non-ideal effects in the equation of state arc partially lifted for Rp > 2Rj. 
Based on the results of a limited number of models computed by BLM for Mp = 0.63 and IMj 
planets around HD 209458 and Ups And ^, we find that 7 crosses the critical value of 5 when 
Rp = R^^ 2Rj in both cases and = ^-c = (10~^, lO^^-^L©) for Mp = (0.63, IMj) respectively. 



Although 7 continues to decline for Rp > Rc, we introduce a power-law approximation 

[ Rp \ ^ 



Co \2R 



(41) 



for a limited range of Rp. In the range Rp = 2 — 5Rj, we fit the results of BLM with a 
normalization coefficient Co = 1O~^L0 and 7 = 3 for a Mp = 0.63Mj planet. This prescription is 
chosen for computational convenience and consistency, i.e. C^ = Cc for Rp = R^'^ 2Rj. In §5, we 
present additional numerical models which also indicate 7 < 5 for i?p > 2Rj for other values of 
Mp. 

In this limit, the thermal equilibrium is linearly unstable to tidal dissipation. A small 
enlargement in the planet's size from R^ would cause tidal dissipation to increase more rapidly 
than the radiative losses on its surface. However, e decays in the meantime. The question of 
whether the planet's Rp would finally reach the Roche radius depends on the magnitude of e: there 
exists a threshold e beyond which the planet is able to be infiated to its Roche radius, meaning 
that this process relies on a nonlinear instability. Similarly, a small reduction in Rp from Re would 
lead to rapid cooling, and the planet would contract until Rp is reduced to Rc. The characteristic 



^The semi-major axis of HD209458b is actually about 0.045 AU. The only difference modeled by BLM for giant 
planets located at different a is the degree of stellar irradiation. However, BLM show that inflation of giant planets 
due to stellar irradiation is negligible compared with tidal inflation. Therefore the estimate performed in the text 
section based upon replacing 0.045 AU by 0.04 AU should serve as a good approximation. 
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contraction time scale at Rp = Rc is 



3.3. Necessary and sufficient condition for runaway planetary inflation 

The above criterion for tidal inflation instability (7 < 5) is applied to a state of thermal 
equilibrium. In general, the planets may not be in such an equilibrium, and during the evolution 
toward it, tidal dissipation may modify their orbital properties. Therefore, it is essential to identify 
1) the sufficient kinematic properties which would enable planets to inflate before their eccentricity 
is damped out and 2) the necessary initial conditions which may lead to an equilibrium radius 
Re which exceeds a critical radius Rc{^ "^Rj) for the onset of tidal inflation instability. In the 
determination of the stability criteria, we must take into account the dynamical evolution of the 
planets' orbits. 

Prior to reaching an equilibrium value of i?e, the planet's radius Rp adjusts at a rate in 
accordance with equation (35). Since -Egpin < -Ebindi this equation implies that Rp inflates on a 
time scale 

S,^_*Gi!5 (43) 

-"■p -'-'tide-'''p 

in the limit that the tidal heating greatly exceeds the radiative cooling. 

Note that equation (35) is based on the implicit assumption that Rp is uniquely determined 
by a given set of Rp, Mp, and E^de regardless of the thermal history of the planet's interior. 
In principle, the evolution of stratification of planets due to thermal imbalance in a state of 
quasi-hydrostatic equilibrium should be a function of the heating mechanism. But, gas inside 
Jupiter-mass giant planets is partially degenerate and non-ideal such that the total pressure is 
insensitive to temperature changes. BLM presented several models with different dissipation 
prescriptions. The results of all their models show that a large fraction of the energy dissipated 
within the planet leads to temperature increases without significantly modifying the pressure 
and thermal stratification of the convective interior of the planet. Only the tenuous radiative 
zone in the outermost 10% of the planet's radius responds to heating like an ideal gas. Thus, 
tidal dissipation only affects the boundary condition at the interface between the convective and 
radiative zones as well as the structure of the radiative zone, and our approximation is justified. 

During the process of synchronization, although the total amount of energy available {Espin) 
is limited, it is dissipated on a rapid time scale of rn such that Ende — Em- From equations (28) 
and (43), we find 
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where Jljtj, is the breakup spin frequency of the planet, and trq is the planetary-inflation time 
scale in the limit where only the effect of synchronization is considered. Since ftp < i^^p and ap 
is relatively small, trq > tq and the planet does not have time to inflate in size before its spin 
become synchronized with its mean motion. 

Energy is also dissipated within the planet during e damping. In order for a planet to be 
inflated by e damping against its own gravitational binding energy, e should decay on a longer time 
scale (Te) than the thermal time scale associated with inflation (tr); that is, the system should not 
run out of its dissipation energy source before it overflows its Roche radius. Prom equations(29) 
and (43), we find that, in the limit that tidal dissipation mainly leads to circularization, 

(45) 

The necessary condition for planetary inflation (tr < Tg) is satisfied in the limit that its e is 
greater than a critical value 

e > ~ f^^^V^' = 18 f _Jh^MLM^^_^^\ (46) 
\ (5 M,Rp) 1^0.75/3 Mj Rp 0.04AU ) " ^ ^ 

In the above expression, decreases with Rp so that the most stringent requirement for a 
protoplanet to inflate is set prior to its expansion. If short-period planets migrated to their present 
location within a few Myr after their dynamical accretion of ambient gas in the protostellar disk 
has ceased, their initial Rp would be ~ 2Rj (Pollack et al. 1996). In equation (46), we use a 
fiducial value of 2Rj for Rp which is also comparable to the value of Rc- During the same early 
epoch, main sequence G dwarf stars are rapid rotators and have relatively small Q'^ so that for 
relatively large Rp, (3 is slightly less than unity. Note that in equation (46) cr is independent of 
Q'p such that it is not subject to the uncertainties in the tidal dissipation rate. 

Short-period planets with initial e > cr would inflate until they reached a thermal equilibrium 

with a radius Rp = R^. Those planets with e » cr would retain their initial e as they approached 
the thermal equilibrium. But those planets with an initial e marginally larger than er would 
undergo a substantial e damping prior to reaching the thermal equilibrium. Thus, it is possible for 
some planets to initially satisfy the necessary condition for tidal inflation, and become stabilized 
after their e is substantially reduced. Thus, the instability criterion in equation (46) should be 
regarded in a qualitative sense. 

Equation (37) indicates that Rg is a function of Ende- Assuming the planet has already 
attained spin-orbit synchronization (so that E^ide ~ -Ste), we find from equations (18), (41), and 
(37) that 

Re (eL{Rs)\^ 



R 



■s 



(47) 
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where Rg is a fiducial length scale and 

1/2 



aTe{Re = Rs)Ce{Re = R. 

eL{Rs) = 



GM^Mp 



where we adopt the approximation in equation (41) for C and scale Rg with 2Rj. Note that 
although Mp does not appear directly in equation (48), both cl and Co are slowly increasing 
functions of Mp. 

In the previous section, we showed that a planet would adjust to a stable equilibrium size if 
Re < -Rc(~ "^Rj)- But, for i?e > -Rc) 7 < 5, and the planet is unstable such that it would undergo 
runaway inflation if its current size Rp < Rg or contraction to Rg if Rp > R^. In principle, we 
should set Rg = Rp in equations (47) and (48). But, the critical should not be a decreasing 
function of Rg as indicated in equation (48). For example, the Roche radius of a 0.63Mj planet, 
located at a = 0.04AU around a 1.1 Mq star (HD 209458), is Rl = 4.92Rj. An extrapolation 
from Fig. 3 in BLM shows that an internal heating rate of at least ~ 10~*^'^ Lq is required for a 
coreless Mp = 0.63Mj planet to attain Rg = Rl in this case. Nevertheless, the critical eccentricity 
cl for Rg to reach R^ is less than half that needed for R^ to reach ~ 2Rj (i.e. < O.bei) because 
the efficiency of tidal dissipation increases rapidly with Rp. Thus, analogous to the planet's Rp 
dependence in br, the most stringent requirement for a protoplanet to inflate is set prior to its 
expansion rather than the final post-inflation value of Rl. In the limit that the planet's initial 
radius Ri < Rg, the appropriate sufficient condition for the onset of runaway inflation is 

e > CLiRg = Rc)- (49) 

Once again, note that e may be damp during the course of the planet's size adjustment, especially 
for those systems with e{Ri) marginally larger than eji{Ri). In these cases, the value of the e in 
the above expression should be e{Rc) (see further discussions in §6). 

The power-law approximation in equation (41) and its extrapolation in equation (48) are only 
appropriate for Rp in the range of 1 — 3i?c where 7 < 5. In the limit Rg = Rc, j = 5, = Cc, and 



~ 0.06 
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3.4. The onset of planets' inflation to their Roche lobes 

If a planet migrated to the vicinity of its host star with an initial radius Ri < Rc and with 
e larger than both cl in equation (50) and cr in equation (46), it would undergo runaway tidal 
inflation. For a > qhl where 

'qp/13 10^ Mp 10-6^0^^^'^ / NU / 2Rj\^ 



< so that the condition for planetary inflation is e > e^. But, a planet with a < aRL, 
can only inflate if e > cr. For a Mp = 0.63Mj planet (7 = 8 and Co = W~^Lq), url increases 
with Rp. Although it is possible for eL{Rc) > ^r{Rc) and CLiRh) < &r{Rl) at some locations, 
the runaway tidal inflation following the onset of the thermal instability always leads the planet 
to overflow its Roche lobe because both e^iRc) > ^LiRh) and eR{Rc) > eR{Ri) such that the 
sufficient and necessary conditions for runaway inflation become more easily satisfied during the 
planet's expansion. Note that the magnitude of grl is very weakly dependent on that of Q'^. 
Two orders of magnitude uncertainty in Q'^ would correspond to a factor of two difference in aRL- 
Thus, our results are quite robust despite the uncertainties in the magnitude of Q'p. 

Equation (46) indicates that br oc o^/^ and it reduces to ^ 0.06 near the surface of a star. 
However, even in the limit that e is larger than both cr and e^, the planet's Rp is limited by Rl 
which is smaller than Rc if 

0.0133 f^liV'V^-V'V^lAU. (52) 



Mp) \MqJ \2RjJ 

Inside ac, planets cannot expand into the unstable range. Nevertheless, Rg > Rc (and therefore 
Rl) for e > eL{Rc) which is very small at such close range (see eq. 48). In this limit, planets with 
e > (so their e- value may be maintained ) would overflow their Roche radius for all values of 

The metallicity of many planet-bearing stars is higher than the solar value. It has been 
suggested that the outer convcctive zone of some of these stars may be contaminated by the 
consumption of one or more planets which may have undergone orbital decay due to the star-planet 
tidal interaction (Lin 1997). Previous models for such a process focus on the disintegration of 
planets inside the envelope of their host stars (Sandquist et al. 1998, 2002). The present results 
suggest that with a modest e-value, short-period planets may be disrupted during their tidal 
orbital decay prior to entering the envelope of their host stars. 

A planet with a > qrl and Ri < Rc would undergo runaway inflation if e > eiiRc)- The 
expression in equation(50) also indicates that oc a^^/^ so that a planet would inflate and fill its 
Roche radius provided its initial a is less than agL where 
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or initial P is less than P^l where 

Over the possible range of Jupiter's Q-value 5 x 10^ < < 2 x 10^, i?c(- 2i?j), and £c(- 1O~^L0 
for the Mp = O.63M0 planet), tidal inflation from Ri < Rc to Roche lobe overflow becomes 
unattainable for a > 0.08 — 0.13 AU or P > 8 — 17 days. For these values of a, Mp, Rc, and /3 ~ 1, 
equations (18) and (43) imply ~ Tg < a few Myrs. 

It is possible, though less likely, that protoplanets may be brought to the vicinity of their stars 
while they are still in the accretion phase so that their initial Ri is larger than R^. In this case, 
the planets are already in a thermally unstable region. These planets would undergo runaway 
inflation if their R^ > R^ and contraction to Rc if their R^ < Ri- A direct comparison between R^ 
and Ri can be made with equations (47) and (48) by setting Rg = Ri such that 



Re _ f eiiRj 
Ri~\ e 



2 



(55) 



R \ 

eURi) ={-f^] eURc) < eL{Rc). (56) 



where 

For a > a^L, runaway inflation from Ri{> Rc) would occur if e > eiiRi). 

In the above derivation of e^, we explicitly assumed that the planet's spin is nearly 
synchronized with its orbit. Equation (27) clearly shows that a significant departure from the 
ri = n synchronization state would lead to an additional source of tidal dissipation of energy. During 
its inflation, the planet's spin frequency would decrease with a rate tlpc = —2QpRp/Rp = —2^p/tr 
if its spin angular momentum is conserved. Nevertheless, quasi-synchronization would be 
maintained if (ipc is balanced by the tidally induced Clp in equation (9) , i. e. 

TO. ~ tr/2. (57) 

Under the assumption that most of the energy is dissipated through circularization (an assumption 
to be satisfied self consistently below), equations (19), (45), and (57) imply that the planet's spin 
attains an equilibrium with 



Qp-n _ 14e2 ( RpV 



^Qp V Rl J 

With such a small departure of Q,p from n, energy dissipation associated with the circularization 
process continues to provide the dominant heat source as shown in equation (27). 



3.5. Planets' initial eccentricity 



Extrasolar planets with a ^ 0.1 AU have a nearly uniform e distribution between 0-0.7. All 
of the dozen extrasolar planets with period between 3-7 days and most planets with a ^ 0.1 AU 
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have negligible e's. This dichotomy is probably due to the tidal circularization of short-period 
planets' orbits during the main sequence life span of their host stars. All the observed short-period 
planets have a > 0.04AU. Although some of these planets may have undergone tidal inflation if 
their initial e > (see eq. 46), but outside ~ 0.08 — 0.13 AU, rapidly increases above unity 
(see eq. 48) so that the rate of tidal dissipation is so weak that the planets' thermal equilibrium 
size is always less than their Roche radius. 

Interior to 0.04 AU, however, a planet would be inflated to overflow its Roche radius provided 
Tr < Te OT equivalently e > Cr which is < 0.18 (see eq. 46). We note that if the extrasolar 
planets with a ;^ 0.1 AU retain their ab initio e distribution, most planets, including those with 
a < 0.04AU, would have initial e> Cr- 

In addition to the extrapolation of the e distribution from the long-period planets, there are 
several theoretical reasons to expect the initial e of short-period planets to be much larger than 
their present values. These short-period planets are most likely to have formed at much larger 
distances from their host stars. Long-term dynamical instability may induce the mature planets 
to undergo close encounters with each other, causing some planets to be scattered to the close 
proximity of their host stars with large e's (Rasio & Ford 1996, Ford 1996, Weidenschilling &; 
Marzari 1996, Lin Sz Ida 1997). In comparison with the value of e^, these planets would satisfy the 
necessary criterion for inflation when their a is reduced (along with their e) by the circularization 
process (Rasio et al. 1997) to a magnitude <0.04AU. 

Short-period planets may also have been brought to the proximity of their host stars during 
their formation stage as a consequence of the planet-disk tidal interaction (Lin et al. 1996) 
through wave excitation at the planets' corotation and Lindblad resonances. Energy and angular 
momentum transfers via the planets' corotation resonances damp their e's whereas those through 
their Lindblad resonances excite their e's (Goldreich & Tremaine 1980). Both processes induce 
gap formation in the disk near the planets' orbit. Massive planets with Mp > O.OIM* clear a wide 
gap such that regions containing all the corotation resonances arc void of gas. The interaction 
between the disk gas and these massive planets through the Lindblad resonances excites e to 
several times the aspect ratio of the disk (Artymowicz 1992, Papaloizou et al. 2001, Goldreich &: 
Sari 2002) which is typically > e^. 

In contrast, modest-mass planets with Mp < O.OIM* open narrow gaps, and residual disk gas 
may still occupy their corotation resonances close to their orbits. The dominance of the corotation 
over the Lindblad resonances induces e damping to negligible values. However, short-period 
protoplanets' e's may also be excited by the processes which halted their orbital migration. For 
example, some young stellar objects are observed to have strong magnetic fields (Johns-Krull et al. 
1999) and their magnctospheres may clear an inner cavity in their surrounding protostellar disks 
(Konigl 1991; Shu et al. 1994). As they enter into the cavity and retreat inwards from the disk's 
inner edge, the planets' corotation resonances evolve into the tenuous regions prior to their more 
distant Lindblad resonances. Thus, during the decline of the planet-disk tidal interaction, the 



e-damping process diminishes more rapidly than the e-excitation process. Eventually, the planets' 
inward migration is terminated when all of their lowest-order Lindblad resonances have retreated 
into the cavity. But, during the transition stage, the planets' e's may grow to be comparable to 
the fractional difference between the corotation and Lindblad resonances which is ^ inside 0.04 



Planets' migration may also be stalled if their loss of angular momentum to the disk is 
compensated by a transfer of angular momentum from the spin of their rapidly rotating host stars 
through their mutual tidal interaction. This process is effective in the limit that a is sufficiently 
small (< Ohait) for the stellar tidal force to be intense and that the stellar spin frequency $7* > n 
(see eq. 22). Fast fi* are observed among young stellar objects, with some reaching the stellar 
breakup values (Stassun et al. 1999). In equation (13), the star's contribution to e is neglected. 
However, the stellar tidal effect would excite the planets' e if /3 < or equivalently. 



Another mechanism to stop planets' migration is through the rapid depletion of the disk which 
would quench their tidal interaction. Signs for additional planets are found around more than half 
of all the host stars with known planets. The e's of multiple planets around the same host star are 
modulated as they exchange angular momentum through secular interaction. During the depletion 
of the disk, the procession frequencies of the planets also evolve with the gravitational potential of 
the depleting disk. This evolution leads to a sweeping secular resonance (Ward 1981) which may 
also induce the excitation of eccentricity for some planets (Nagasawa et al. 2002). 



Planets with e > e,. inflate before their orbits are circularized. If their e > e^, the expansion 
of their envelope would continue until it reached the Roche radius. In this section, we consider the 
flow pattern associated with the mass loss process. 



AU. 




(59) 



which ~ 8n in the limit that Q' 
star (DLM). 



'5 = Qp and the average density of the planet equal to that of the 



4. 



Mass Loss from Inflated Planets 



4.1. Density scale height near the Roche lobe 



In a quasi-hydrostatic equilibrium, the planet's surface is determined by the equipotential 
surfaces. For planets with circular orbits, the dimensionless potential in a frame which co-rotates 
with their mean motion is reduced to 
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where ri and r2 are the dimensionless distances from the host star and the planet, respectively 
(Murray & Dermott 1999). The physical values of ri and r2 are obtained by multiplying them by 
a, and that of U can be obtained by multiplying it by GM^/a. 

For this potential, there axe two locations, LI and L2, near the planet which are saddle points. 
Both points lie on the line joining the planet and the host star, with LI between them and L2 
farther away from the host star than the planet. The distances of the LI and L2 points (-Dli and 
Dl2) to the center of the planet are 



Dli,L2 ^ 



a2 



(61) 



where a = {Mp/3M^y/^ . Gas that expands beyond the LI point falls toward and is accreted 
by the host star, whereas that which expands beyond the L2 point spirals outwards to form an 
excretion disk. 

The cquipotcntial surfaces containing LI or L2 points are distorted from spherical symmetry 
about the planet. On the planet's "night side" (in the direction of L2), the distance between the 
planet's center and the closest point (LI') on the equipotential surface which contains the LI point 
is 

L>^i ~ Dli - ^«'/'a. (62) 
The distance between the LI' and the L2 points is 

AD = DL2-D'L,c^'^ia'/'+a') (63) 

which for a Jupiter- mass planet is of the order of a few percent of -Dli- It is possible for the 
planet's envelope to reach the equipotential surface containing the LI point but not that containing 
the L2 point. 

However, it is also possible for both equipotential surfaces to be contained within the planet's 
atmosphere. For an isothermal planetary atmosphere, the density {p) scale height 

h - " =^ (64) 

The gradient of the potential near the LI' point is 

„^ GMtdU , s/oGMt 
ar2 

which is somewhat softer than that of isolated planets. The ratio 

^ 9 TlgTpa ^gg^ 



AD 8 figOMp 

where Tp and Cs{rp) are, respectively, the temperature and sound speed of the planetary surface, 
and TZg and /x^ are the gas constant and molecular weight. On their day side, planets are heated 
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to an equilibrium temperature Tp ~ T^{R^/a)^/^, where the surface temperature of solar-type host 
stars is ~ 6 X 10^ K. For Jupiter-mass planets at ~ 0.04 AU, hp ~ AD, such that both LI and 
L2 may be contained in their atmosphere as the envelope fills the Roche radius. 

In the above estimate, we used the day-side equilibrium temperature to estimate Tp. But 
both LI' and L2 points are on the "night-side" of the planet and are shielded from the stellar 
radiation. In principle, the night side of a mature planet could have a substantially lower surface 
temperature if the cooling time scale there is much shorter than the time scale for sound waves to 
cross the dark hemisphere. But, the resulting longitudinal pressure gradient would excite large 
scale circulation with speeds comparable to Cg of the hot side (Showman &; Guillot 2002, Burkert 
et al. 2002); it is not clear yet what the resulting temperature on the dark side would be. For 
planets with relatively young ages {tp ~ 10^ — lO"^ yr), however, the effective temperature of 
unheated Jupiter-mass planets arising from their intrinsic luminosity (Burrows et al. 1997) falls 
in the range 500 - 800 K. If the planets migrated too close to their host stars and became tidally 
inflated shortly after they were formed, the stellar heating would introduce little temperature 
difference between the day and night sides of the planet. 

In addition, if there is sufficiently intense energy dissipation in some planets to inflate their 
Rp beyond their Rl, their surface luminosity Cc{Rl) would exceed 1Q~^Lq (see §3.2). For these 
tidally inflated planets with Rp ^ Rl > 2i?j, equations (40) and (41) are the more appropriate 
approximations. Using a black body scaling law C/ Lq = {Rp/ Rq)'^{Tp/Tq)'^ where Tq = 5800K is 
the effective temperature of the Sun and equation (41), we find that when Rp = Rl, 

which for a Mp = 0.63Mj planet is ~ 700 — 1000 K. In this limit, hp may be a fraction of AD. 
But, we also note that near both the LI and L2 points, the magnitude of V* vanishes so that hp 
may be substantially larger than that in equation (66) for the LI' point. 



4.2. Mass flow induced by envelope expansion 

We consider the flow across equipotential surfaces below but near the LI point. For 
computational convenience, we adopt a spherically symmetric approximation for the equipotential 
surface. During the quasi-static expansion of the envelope, the mass flux advected from the 
planet's interior across U is 

Ma{U) ^4TTp{U)r^{U)Vr{U) (68) 

where Vr{U) is the outflow speed across the equipotential surface located at a radius r{U) from 
the planet's center. Vr should not be confused with the size expansion rate Rp. Assuming the 
planet's expansion is homogeneous, we can approximate 

mn^A^A. (69) 
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where A is a positive dimensionless structure parameter. We show below that for a homogeneous 
interior, A = 3 and it is of the same order but shghtly less than 3 for more realistic models. 

Neglecting departure from spherical symmetry, we define the equipotential surface U to be 

located at a distance Riit) from the planet's center. At any given time ti, the mass interior to 
Rl is Msiti) = p{r,ti)r'^dr. For Rp > Rl, the mass outside i?^, i.e. that between Rl 

and Rp, would be lost. At a later time t2 = ti + At, the planet's density distribution, size, and 
Roche radius evolve at rates p, Rp, and Rl respectively. The difference in the mass contained in 
Rl becomes 

M.(,,) = ^ = ^MlM_mA^ ^ ,MR.(h))Rdt.n, + A. t^ir. (70) 

In a Roche potential, 

Rl 1 / mh) \ 

Rl ^W!^-^''^prHr)' ^ ^ 

Prom equations (70) and (71), we find 

In the limit that Rp is slightly larger than Rl^ mass loss occurs near the atmosphere where the 
density is low and the second term on the left hand side of equation (72) is very small. But, after 
the planet's atmosphere is removed, the envelope continues to expand. The atmosphere merely 
responds with an expansion to restore a state of hydrostatic equilibrium. The actual value for 
p{Ri,ti) is that prior to the removal of the envelope between Rl and Rp. 

An approximate solution can be constructed in which wc consider the limit that Rl is much 
smaller than the original Rp. We first consider the rate of change of density in terms of the rate 
of change of some intermediate reference radius Rr- The results of numerical models of planetary 
interiors (see §5) indicate that the inflation is nearly homologous, i.e. the ratio of the half density 
radius changes as R1/2 oc R^^ where ("1 ~ 0.47. The central density pc oc R^^ where ("2 — —1.35. 
These numerical results also indicate that the density within R1/2 is nearly homogeneous such 
that the mass contained within it ~ 47TpcRy2/'^ approximately constant. 

In order to derive the rate of change in density p, we first consider a reference radius R^ 
which, at some time ti, contains a mass 

/o^"^*'^ p{ti)r^dr = Mr. At t2 = h + At, the radius which 
contains the same amount of mass has evolved to Rr + i?r At while the density evolves at a rate p 
so that + p{ti)At)r^dr = Mr. Note that although Mp may decrease due to 

mass loss, but, Mr is chosen to attain the same value. Prom these expressions of Mr, we find 



pihydr = -p{Rr,t,)RrihfP^^. (73) 

Rr{ti) 
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Note that the rate of change Rr includes that due to tidal dissipation as well as mass loss. For 
isolated planets with Mp > Mj, dRp/dMp < because the planet's interior is in a partially 
degenerate state. Thus, mass loss may also lead to size increases. 

Without the loss of generality, we can choose Rr to he Rl in the above expression so that we 
find from equations(72) and (73) that 

Mp \ SMp ) Mj, ' " ' ■ ^ ' 

Setting Ms = Ma, we find from eqs (74) and (69) that 

_ 47rp{RL)Rl {Rr/Rr)R^ 

Mp Rp/Rp 

Shortly after overflowing the Roche radius Rl, {Rr/Rr)Ri^ — Rp/Rp and p{Rl) is very small. But 
in the advanced stage of Roche lobe overflow, {Rr/Rr)Rj^ — R1/2/R1/1 — CiRp/Rp but p{Rl) 
may be only slightly smaller than pc- For a nearly homogeneous planetary interior A ~ 3 and if 
it is more central condensed, A < 3. Our numerical models for their interior indicate that under 
the heating due to the tidal dissipation, planets establish a quasi hydrostatic equilibrium prior to 
reaching a thermal equilibrium. Since most of the gaseous giant planets' interior is convective, 
the temperature gradient is essentially adiabatic and the average global adiabatic index is ~ 1. In 
these models, A is of the same order but slightly less than 3. The numerical results to be presented 
in §5 do not depend sensitively on the actual value of A. 



4.3. Mass flux through the LI &; L2 points 

As the planet's envelope expands and fills its Roche lobe, gas in its atmosphere encounters 
and flows through the saddle-shape nozzle in the equipotential surface. Near this region, the above 
approximation based on spherical symmetry can no longer be applied. Instead, we consider a 
cylindrical geometry with an axis of symmetry along the line joining the planet and its host star. 
The planet's mass loss rate through the LI point between the equipotential surface which contains 
it and that with an arbitrary value Un is 

MLiiUn) ^ / p{U)csiU)—dU (76) 
JUi dU 

where a is the cross section of the nozzle (Lubow & Shu 1975, Pringle 1985). Note that the total 
mass flux can be obtained by setting Un = Uoo- In a coordinate system in which x is defined to 
be the direction joining the centers of the planet and its host star, y is orthogonal to x and in the 
orbital plane, and z is orthogonal to both x and y, the narrowest opening of the nozzle occurs at 
X ~ xli, with a = 7r(sa)^ lies in the y — z plane perpendicular to the x axis and sa is the radius 
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of the nozzle. The dimensionless variable s is a function of U and it vanishes for U{r < Dli). In 
the neighborhood of LI, dU/ds 2± 3s so that 

da 27ro^s 27ra^ 



dU {dU/ds) 3 
For the same region, we find from equation (64) that 



(77) 



/ 3GMts2\ 



where pi is the density at the LI surface, and from equation (76), we find 



where Sn = [2(f/„ - C/i)/3]V2. 



1 — exp 



2acl 



(79) 



The expansion of the planet's envelope would continue so long as Ma{U) > Mli- As the 
planet's envelope climbs out of the potential well, p would increase which leads to enhanced 
spillage via the LI point. Wc now estimate the critical condition which may lead to Roche lobe 
overflow via the L2 point. The difference in the potential between LI and L2 points is 

MJ ~ 2a^ (80) 

so that SL2 = 2(Mp/M*)V2/3 and 



Mli{sl2) 



27rpi a'^c; 



3^3 r 



1 — exp 



2GMp 
Sac^ 



(81) 



3GMt 

The expansion of the envelope would terminate if a mass transfer equilibrium can be established 
in which 

Mli = -Ma{U{Ll)) (82) 

if we neglect the adjustment of the equipotential surface due to mass loss and orbital evolution 
(for further discussions, see next subsection). 

Prom equations (69), (81), and (82), we find that the mass transfer equilibrium via the LI 
point would only attainable if the local density 

where the planet's average density within its Roche lobe Pa = 3Mp/47ri?| ~ 10~^g cm~^ for a 
Jupiter mass planet at 0.04 AU. Since the envelope of the inflated planet expands at a speed 
Rp « Cs, p{Ll) « Pa- For Rp = Rl and e ~ 0.14, equations(43) and (83) lead to 

7A/3e2 fGMp\ f RLn\ ^ . ^ o , 
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Since pa oc a ^, Cg oc a and n oc a Pc is proportional to a ^^/^e^ and it would exceed 
10~"^g cm~^ at a < 0.025AU, which corresponds to a 1.5 day period around a 1 M© G dwarf. 

When planets overfill their Roche radius, p{Ll) is comparable to the density patm at their 
efi^ective radius where the optical depth is of the order unity. Although V'l' = at the LI point, it 
is ^ GMp/R\ or that given in equation (65) elsewhere on the equipotential surface which contains 
that point. At a = 0.04 AU where the day-side temperature of 3], ~ 1 — 2 x lO^K, we find from 
equation (64) that /ip ~ 2 — 3 x 10^ cm for a Jupiter-mass planet. For this Tp, the opacity of 
the atmospheric gas is k ~ 10~^ cm^ g~^ so that patm ~ 10~^g cm~^. For modest e values and 
a > 0.04 AU, patm > Pc such that the optically thin atmosphere engulfs both the LI and L2 points. 
In this limit, the density at the L2 point can be derived from equations (64) and (80) to be 

On a Mp = IMj planet at 0.04 AU from a G dwarf where Tg ~ lO^K, p2 ~ 0-2pi so that gas near 
the Roche lobe mostly fiows toward the host star via the LI point. 



From equation (84) we find Pc > Patm 

for planets with 

a < aL2 



/ Sp, \ rGMp\ / 7X/3 ^ fR,7u ^ ^2 
\PatmJ \clR„J \6g„Q', 



4/15 

R* (86) 



where p^ = 3M^/AttRI, = (GMjRl)^/^, and = {UgT^pgY'"^. For e ~ Ci? ~ 0.14, 
ol2 — O.OSAC/. Equivalently, with a given a, pc > Patm for planets with 



e > eL2 




1/2 

(87) 



At a = 0.03 AU, eL2 — 0.14. In this limit, the opaque regions of the envelope extend above the 
equipotential surface which contains the LI point. The temperature and density scale height in 
the envelope are greater than their surface values so that the density differential between the LI 
and L2 points is reduced. Consequently, gas at the planets' effective radius may flow through both 
LI or L2 points with comparable rates. An accurate determination of the ratio of the mass flux 
and the specific angular momentum fiux through the LI and L2 points needs to be analyzed with 
more detailed numerical hydrodynamic simulations (Burkert &; Lin 2003). 

In §3, we showed that the necessary and sufficient criteria for planetary inflation are e > 
and e > cl- The instantaneous tidal radius of eccentric planets may modulate by a fraction e. At 
the periastron, gas outside the planet's minimum tidal radius is removed. During other phases, 
the planet's atmosphere is refilled as it adjusts, on its internal dynamical time scale, to an evolving 
quasi hydrostatic equilibrium, while the underlying envelope re-expands on a much longer internal 
heating time scale. If the ratio 

hp ^ (3a) -ji^Tpa 
2eDLi ~ 4e p.gGMp 
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is comparable to or larger than unity, the planet's atmosphere would overflow the Roche lobe at 
all orbital phases and the mass loss would be continuous. Otherwise, bursts of mass loss would 
occur during each periastron passage. 



4.4. The rate of planets' mass overflow through their Roche lobe 

Planets' Roche lobe overflow also leads to their orbital migration. If the overflow occurs 
through both LI and L2 points, angular momentum exchange between the infalling and outflowing 
transferred gas would result in orbital migration. In general, mass loss would lead to a reduction 
in the Roche radius at a rate 

Rl^( (Mli+Ml2) a _ 

Rl \ SMp a 1-e) ^ ' 

along with the equipotential surfaces which contain the LI and L2 points. The shrinkage of the 
Roche lobe would enhance mass flux across the equipotential surfaces which contain both the LI 
and L2 points by an amount 

such that the effective mass flux across the equipotential surface of the newly adjusted Roche lobe 
U'{Ll) is 

Me(C/'(Ll)) = Ma(C/(Ll)) + Mr. (91) 
In a mass transfer equilibrium, this mass flux is lost through the LI and L2 points such that 

(Mil + Ml2) = -Me{U'{Ll)). (92) 

For e > eL2, mass is lost via both the LI and L2 points. In this limit with the simplified 
assumption d ~ ^, we find, from equations(69), (89), (90), (91), and (92), the planets' mass loss 
rate to be 

For planets with a relatively flat density distribution (i.e. A is smaller but close to 3), the 
magnitude of Mn can greatly exceed that of Ma{U{Ll)). 

The reduction of the planet's mass reduces the radial velocity amplitude K for Earth-based 
observers which is determined by 

^ ^ MpSinz j^ 2^G(M,) y/^^^ 



^This over-simplified assumption is just for the demonstration of this possible scenario, and will be used in the 
next section. What exactly happens to a in this limit is undetermined in this paper as stated in the preceding section 
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where P is the planet's period and i is the inclination angle of the orbit with respect to the line 
of sight. For an observational limit of 1 m s~^, planets with a few Earth masses and -P ^ 3 days 
would be undetectable with the Dopplcr method. Thus, a substantial mass loss may account for 
the absence of any extrasolar planets with a <^ 0.04AU. 

If, however, e < 6^2 , planets would overflow their Roche lobes mostly via the LI point and 
they would move outwards as they absorb the angular momentum and energy from the transferred 
gas which flows toward their host stars. This dynamical property is analogous to conservative 
mass transfer from a low-mass to a high-mass star in an interacting binary system. Assuming 
all the transferred material is accreted by the host stars and the total angular momentum of the 
system is conserved, the semi-major axes of these planets would expand at a rate 

^ = 2-^-2^. (95) 

Unless the planets have very small a's, the matter lost from their Roche lobe cannot directly 
strike the surface of the host star, and it forms a circumstellar disk. During its infall onto and 
subsequent viscous evolution within the disk, both energy and angular momentum are tidally 
transferred between the overflow gas and the planets' orbits. Unless the net energy and angular 
momentum transfer rates are exactly balanced, eccentricity of the planet may be modified at a 
rate e^/e = r/M^i/Mp. The actual magnitude and sign of rj is determined by the combined and 
competing effects of corotation versus Lindblad resonances (Goldreich k, Tremaine 1980). Despite 
the mass transfer, the binary stars' orbits in cataclysmic variables remain circular, which implies 
that the magnitude of is negligibly small compared with that due to the tidal circularization 
in equation (13). For the remaining discussions, we neglect the contribution associated with e^. 

The semi-major axis of an eccentric orbit would shrink due to tidal orbital circularization 
at a rate Oe = 2aee/(l — e^) (see eq. 21). This effect of circularization on a is more prominent 
for larger e, which would be the case for extrasolar planets. The transfer of angular momentum 
between the planet's orbit and the stellar spin also introduces = laJ^jJ^ (see eq. 20) which 
is generally much smaller than that induced by mass loss or orbital circularization. Finally, the 
planets' tidal interaction with their residual nascent disk also contributes to inward migration on 
a time scale t^- Consequently, the total change of a comes from these three effects: 

^ = ^^ = -2^ + 2-^ + 2:^-1. (96) 
a a Mp (l-e2) Jo 

The above equation, together with equation (89), implies that the Roche lobe changes at a rate 

Rl 5Mli e 



Rl 3Mp 1 + e 



(97) 
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where we have neglected the effects of stellar spin synchronization and planet-disk tidal interaction. 
The first term on the right hand side of the above rate equation indicates that the retreat of 
the planet's orbit due to mass overflow would enlarge the Roche lobe. On the other hand, the 
shrinkage of the planet's orbit due to orbital circularization would reduce the Roche lobe in size 
as is indicated by the second term. 

In a mass transfer equilibrium, we find, from equations (69), (90), (91), (92), and (97), the 
mass transfer rate via the LI point to be 

Since A is positive and greater than unity, the magnitude of Mli is significantly reduced from that 
of M„(C/(L1)). 



5. Mass and dynamical evolution due to Roche lobe overflow 

The discussions in the two previous sections indicate that the evolution of Mp, e, and a are 
closely connected to each other. In this section, we present analytic approximations to assess the 
dynamical evolution and the outcome of planets which undergo tidal inflation instability. 



5.1. Five stages of planets' tidal inflation 

For illustration purposes, we divide the evolution into five stages. 

Stage 1: Runaway inflation. Wc consider a planet with an initial semi major axis at, 
eccentricity Cj, and radius Ri, so that its initial angular momentum is Jj = Mp(GM^{l — e?)ai)^/^. 
These parameters are determined by processes which brought the planets to the vicinity of their 
host stars. For example, during the epoch of planet formation, young solar type stars may have a 
radius 2-3 Rq so that the planet-star and planet-disk interaction would lead to a halting radius 
in equation (22) = i?hait ~ 0.01 — 0.05 AU. Or if planets are scattered to the vicinity of their 
host stars, their and ej may be relatively large but their periastron distance may be only a few 
stellar radii. 

Provided i?j < i?^ ~ 2i?j, the necessary and sufficient conditions for tidal inflation instability 
are, > eR{Rc) > eiiRc) for < aRL{Rc) (see eq. 51) and > eL{Rc) > eR{Rc) for Oj > url- 
For Ri > Rc, these conditions are > eR{Ri) > eiiRi) for < aRL{Ri) and Cj > eL{Ri) > eR{Ri) 
for Oj > aRL^Ri)- Unstable planets with ai < Qc overflow their Roche lobe prior to runaway 
inflation whereas those with Oj > Oc undergo runaway inflation to overflow their Roche lobe. 

These criteria for runaway tidal inflation are derived under the assumption that e is not 
damped significantly as the planets adjust their sizes toward Rc- This approximation is generally 
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adequate for e >> e^. But if e is marginally larger than initially, it may be damped by an 
non negligible amount and become stabilized before Ri can reach Rc. Such a possibility is best 
analyzed numerically (see §6). 

The dimensionless expansion rate 

is obtained from equation (35) by interpreting Rp as the adjustment rate of the planets' half-mass 
rather than the Roche radius. Note that strictly speaking, in the above equation should not 
be the value related to Ce as shown in eq(48) but the value related to £.. However, our internal 
structure code indicates that C o::^ Ce for a given Rp and Mp, and we adopt this assumption 
hereafter in this section. Equation (99) is appropriate for all Rp and a's. For e > er, the planets' 
radius adjustment time scale, r/j = Rp/\Rp\ < Tg at the onset of tidal inflation. Since r/j is a 
rapidly decreasing function of Rp and the planets' semi major axes evolve only after they overflow 
their Roche radius, the duration of stage 1 is determined by tr with Rp = Ri such that 

. r.. Ri f 4 \ ( Q'p\ f Mp\ f Mq\^I^ f a V^'^ f2Rj\\^ 

- - ^^-""^ - RM^) - [^) [w] ht) hi) (^) (i^) 

(100) 

which, for Ri Rc ^ ^Rj, is longer than the expected values of Td{'^ 10^~^ yrs). Thus, it would 
be a reasonable approximation to neglect the effect of planets' tidal interaction with the disk if 
they have already terminated their accretion and contracted to < 2Rj when they migrated to the 
proximity of their host stars. However, if the protoplanets' orbital migration, induced by their 
interaction with the disk occurs on a time scale which is comparable to or shorter than that for 
them to to contract (Trilling et al. 1998), they may reach the vicinity of their host stars with 
Rp ~ 2Rj. 

Equation (100) indicates that stage 1 could be completed for those planets which migrated to 
a < 0.04 AU while their host stars evolved through the classical and weak-line T Tauri phases. At 
the end of stage 1, the tidally inflated planets lose mass through Roche lobe overflow (see stage 
2 below). All the metal-rich planetary tidal debris flowing out of the LI point is accreted by the 
host stars and thoroughly homogenized in their surface convective envelopes. Early planetary 
disruption and stellar accretion is unlikely to significantly modify the metallicity of the host stars' 
convection zone because it is extensive. But ri > 10^ yrs for planets which terminated their 
migration at slightly larger radii (at a > 0.05 — 0.06AU). In this limit, the onset of the planets' 
Roche lobe overflow is delayed until their host stars have fully evolved onto the main sequence. 
For example, the total mass of the surface convection zone of a G dwarf host star is reduced to an 
asymptotic value of ~ O.O2M0 after ~ 30 Myr (Ford et al. 1999). Thereafter, the accretion of 
the metal-rich tidal debris from the Roche-lobe overflowing Jovian-mass planets may significantly 
contaminate the outer regions of their host stars. Most planet-bearing stars are observed to have 
metallicity higher than the field stars (Gonzales &; Laws 2000). 



(99) 
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Stage 2: Orbital evolution of Roche-lobe filling planets (e^ — e\{RL) > e|j(i?L)). After filling 
its Roche lobe, an inflated planet continues to expand and to lose mass through both LI and L2 
points at a rate 

^ = ^f^ + ^V (101) 

Mp (l-A/3) \Rp 1-eJ ^ ' 

with d ~ if their Qi < aL2 (see eq. 93). But, for ai > aL2 the mass loss and migration rates are 

MLi = ^ (^ + ^] (102) 

a ee Mr-i 

- = ^ - 2^ 103 

a (1 - e2) Mp ^ ' 

(see eqs. 98 &: 96). Note that negative values of M^i correspond to outflow. Since accretion 
cannot occur in this context, Mli must be set to zero if the right hand side of equation (101) or 
(102) becomes positive. 

Neglecting the effect of e excitation due to planets' interaction with the disk and the 
dissipation within the stellar interior, the rate of e damping is 



e 



2l7r /M„\2/^ 1 



e 32/32g;, VM 



(104) 



for both limits (see eq. 13). The dimensionless expansion rate in equation (99) implies that, 
regardless of the mass loss, the expansion of the inflated planets continues provided their R^. > Rj^ 
(negative values of R correspond to contraction.) Note that for e > e^, the planets' radius 
adjustment time scale, r^j = Rpl\Rp\ < Tg at the onset of tidal inflation. Since eji oc Rp these 
inequalities are preserved as the planets inflate and fill their Roche radius. For planets with 
Rp = Rl, 



which is roughly oc a~^^^ according to eq(103) and therefore is weakly dependent on a. On the 
other hand, 



rp-R^ - 0-056 / Q'p 

eL{ns itL) 26(5-7)/6 I 106 iQ-SLp 



1/2 



/ 

1/2 



Mq ( a Y 

vo-04Auy 



(5+27)/12 



^ 27^ (Q^^^\ - (^^^-)/^ 

~ 263.9 \1Q^10-^Lq) Udays; ^ ^ 

where the rough approximation is made that the planet is in thermal equilibrium. For a 
Mp = 0.63Mj planet (with 7 = 3, £<, = W^Lq, and Rp = Rl), 




eL{a,RL) = 0.02[^\ — — . (107) 
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which is cx a^^^^. 

During stage 2, — e\{RL) > e|j so that Rp and e have opposite signs, i.e. eccentricity 
damping leads to further inflation. Consequently, the planets are inflated to fill their Roche radius. 

In addition, the magnitude of tji is smaller than that of Tg, i.e. the planets' size adjusts faster 
than their eccentricity changes. For large Cj (>> eji{ai) and >> eilai)), equation (99) implies 
negligible changes in e from Cj while a sizable fraction of Mp is lost and a is increased significantly. 
But for marginally larger than eR^Ui) and eL^ai), the planets' tidal inflation and mass loss is 
accompanied by an e reduction by a comparable fraction. 

Note that if mass loss does not significantly modify a, as is the case for < aL2, would 
remain constant. Consequently, the stage of runaway inflation continues until the planets are 
disrupted locally. But for Oj > aL2, most of the mass is lost through the LI point while both a 
and eL^Rh) increase. 

At the end of stage 2, e^(oi,) = e\{RL) + e\{aL,RL) which occurs at a = ol (see expression 
for ol and the corresponding period Pl in the discussion for stage 3 below). While mass transfer 
is sustained, equations (99), (102), and (103) implies that 



a 
a 




Ml a 




(108) 
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Mp (1 + 5A/3) \^ e| 1 + eJ e 

Since A ~ 3, tm and Tq are < tr, but comparable to Tg which is an increasing function of a and 
P. Thus, the duration of stage 2 is 

which is much shorter than the duration of stage 1, i.e. ri c:i Tii{ai,Ri). 

Thus, runaway tidally inflation requires a substantial lead-in time which is followed by a brief 
expansion and orbital migration phase. Because they arc substantially inflated, the theoretical 
probability that the transit of Roche-lobe overflowing planets in front of their host stars are 
observable by distant observers is 

which is larger than for the present day short-period planets due to the enlarged sizes of both 
i?* (for Ti < 10 Myr) and Rp. The actual probability of observing such events may be much 
smaller because 1) only a small fraction of stars have short-period planets and 2) the duration 
of this Roche-lobe overflowing stage 2 is brief compared with the life span of the stars. Note 
that during the transit, the host stars' magnitude is reduced by Amj ^ 51ogio(-Rp/-R*) which for 
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the Roche-lobe overflowing planets may be substantially larger than the present-day short-period 
planets. 

Stage 3: Stagnation and eccentricity damping (e|j(i?L) > — e|,(i?L) > 0). The discussion 
below on stages 3-5 is relevant only for the cases with > aL2- Although the magnitude of the e 
reduction depends on the initial conditions, e|(a, Rl) increases monotonically to — e|j(i?L) such 
that 

ej ( e\aL) - cURl) ] ^ ei(a^,i?^) ^ /a^yV2 / R, ^-^ 



In deriving the above equation, we ignore the slight dependence of Cq on Mp for the simple 
illustration. Note that eR{RL) weakly depends on a (see eq 105) and it is less than eR{Rc) (see 
eq. 46). The onset of tidal inflation instability also requires Cj > eL{ai,Rc) so that the left hand 
side of equation (112) is greater than unity. For a Mp = 0.63Mj planet (7 = 3), 



with a corresponding period 



0.17 (e\aL) - cURl)) ^ ij^) AU (113) 







Pl = 25.5 [e\aL) - e|(i?L))'^'' j days. (114) 

In the limit Cj >> eR^Ri), the planet's e retains its initial magnitude Cj when the above 
conditions are established so that ol ^ 0.17ef ^\m*/M0)1/3aU, and Pl = 25.5e-/^\l06/QP^/^^ 
days. But, for marginally larger than efi(Rc) and cl^Rc), ai, ~ Oj, i.e. the planets would not 
migrate significantly before its tidal inflation is stalled. 

As — e\{RL) — e\{RL) reduces to just below zero, becomes long compared with Te, 
although the planet continues to expand provided e > eL{Ri). But, the reduction of e also leads 
to a reduction in a and e^. The planet would continue to fill its Roche lobe and retreat in a if 
deLiRh) / de > 1. Prom equation (107), we find 



deURL) ^ dcLiRL) da 
de da de 

But, from equations(102), (103), and (99), we find 
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so that with = e|j -I- e| , 

dcLiRL) II ( A 



de 2 U - e2 (1 + 5A/3) (1 + e) 



(117) 
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which, for e < 0.6, is less than unity, primarily due to the planet's retreat in a as a consequence of 
angular momcntTim conservation. Thus, for modest values of e, the quasi-equilibrium cannot be 
maintained and e rapidly declines below e^,. 



During stage 3 when the planet's e decreases from (e|j + e|)^/^ to e^, Rp remains positive but 
tr is larger than Tg. For 

/M^\ 7/33 

Q.lhQ'J \Mj) \Mq, 

eji > cl so that e decreases by a substantial fraction (~ eR/ei)-, Mp decreases and a increases 
by a modest amount. But for a > aRL, cr < cl so that e, Mp, and a do not change significantly 
during stage 3. The duration of this stage is T3 ~ Te{aL, Rl), *-e. comparable to T2. 

Prom equations (102) and (99) we find that Roche lobe overflow is quenched when 

el{RL) -e^ + -^eR{RLf = 0. (119) 
1 + e 

For relatively small eR{RL), the solution for the above algebraic equation is 



e = et ~ CLiRi) 1 + 



eiiRL 



2{1 + eL{RL)) 



eR{Rl 



-I 2N 



cURl) 



(120) 



The termination of Roche-lobe overflow is the result of two competing processes. The reduction in 
e leads to decreases in a but an increase in Rl . For e > e^, the expansion of Rp is faster than R^, 
and equation (103) corresponds to the conservation of total angular momentum. Thus, 

Mp{at,et) = MJ ''f ~i\ ] ' (121) 

which is the asymptotic mass of the planet. For e slightly less than e^, the expansion of R^ is 
faster than Rp. Thereafter, equation (102) is replaced by M^i = and equation (103) is reduced 
to a requirement for subsequent a and e to obey the conservation of specific angular momentum, 
i.e. a(l — e^) = at(l — ef) where at is the semi-major axis when e is reduced to 64. Because the e 
range in stage 3 is limited and has a weak e dependence, ~ o^. 

Stage 4: Onset of contraction {e\{Ri,) > e\{RL) — > 0). When the planet's e decreases 
below eiiRL), Rp and e attain the same sign such that the planet contracts despite the tidal 
dissipation (see eq. 99). Although Rl > Rc so that the Roche-lobe overflowing planets are 
unstable. Re < Rl for e < eLiRi), i-e. the tidal dissipation can no longer supply adequate flux of 
energy to counterbalance the radiative losses. In this limit, the planets undergo contraction. 



The initial contraction proceeds on a time scale \tr\ > Tg. The reduction in e and Rp quickly 
break the thermal equilibrium such that C» — £?tide- In this limit, the contraction time scale 

4 . . 2 

yr (122) 



TeiRp) = [-) re = -^re ^ 1.6 X 10 ^-^ j 



-34- 



for a Mp = 0.63Mj planet (see eq. 42). Equation (102) indicates that the outflow is quenched by 

the contraction of the planets' Rp despite a reduction in their Rl as a consequence of a reduction. 
Setting Mil = in equation (103), the conservation of specific angular momentum a(l — e^) 
implies that a has the same sign as e and that the subsequent circularization leads to orbital decay. 

Since the planet's e is damped before its size is significantly deflated, the duration of stage 
4 is T4 ~ Te^Rh) which is also comparable to both T2 and rs and the age of typical classical and 
weak-line T Tauri stars. For a < grl, eji{a,RL) < eL{a,RL), tr remains to be larger than Tg as e 
diminishes well below eR. In this case, the planet retains its extended -Rp(~ Rl{o,i) 10i2j) until 
its e vanishes on a time scale Te{at,RL) (see eq. 104). The asymptotic value of the semi-major 
axis is a ~ at{l — ef). The planet's extended size and aspect ratio (~ 0.1) makes it detectable 
with transit searches among young stellar objects. For a > url, eL{RL) > CRiRi) and the 
initial contraction evolves into a rapid deflation to Rc when e decreases below {e\ — ej^)^^"^. The 
fractional e reduction is 5e/eL{RL) — &r{Rl)'^ /'^&l{Rl)'^ < 1 with a similar change in a. In this 
case, the duration of stage 4 is ~ (e|j(i?L)/2e|,(i?L)re. 

Stage 5: Runaway deflation and return to stable equilibrium {e\{RL) — > e\{Ri,))- The e 
damping reduces the energy dissipation rate and the planet's equilibrium radius Re to be less than 
Rc or ~ 2Rj. For a < qrl, e vanishes on a time scale TeiRi) « tr{Rl) ~ Tc{Rl)- Thereafter, 
the planet contracts from the Roche radius to Rc on a time scale Tc{Rc)- Although planets located 
at a > qrl may retain a substantial e{r^ gl^Rl)) when their e\{RL) — > e^^^RL), they also 
contract to ~ Rc on the time scale of tr ~ Tc{Rc)- But, since Tg oc Rp^-, the e damping time scale 
for these planets increases rapidly with their contraction. 

Since CR^Rp) oc Rp^^"^ (eq. 46), and eL{Rp) oc R^e ^''^^ (eq. 48), the magnitude of both er 
and cl increase by the contraction of Rp when 7 < 5 or i?p > Rc- As e becomes smaller than cr 
and e^, the effect of tidal heating diminishes monotonically. In both limits, the runaway deflation 
reduces Rp below ~ i?c on a time scale Tc{Rc) ~ 10^ yr which is the duration of stage 5, rs. 

After Rp is reduced to ~ 2i?j, the planet's contraction slows down. Initially, C » —En^e 
so that the primary source of energy for the planet's luminosity is the release of gravitational 
energy. But, this energy release rate decreases with time (Burrows ct al. 1997). If e is damped to 
a negligible magnitude, Rp would contract to become comparable to that of Jupiter today. The 
retention of any residual e would provide an additional energy source such that the contraction 
would continue until an energy equilibrium is established and Re attains the values given by 
setting C = — £^tide in equation (40). 

Associated with this energy dissipation, the planets' e would be dissipated on a time scale 
Te{Re) (see eq. 18). But for relatively large a's and small Rp, Te is likely to be much longer than 
the expected life span of typical solar type stars and e is not signiflcantly damped after stage 4. 
Similarly, the evolution of the planets' a is also likely to terminate. 
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5.2. Semi-major axis expansion and mass reduction factor 

The final dynamical properties of a tidally inflated planet are determined by the conservation 
of total angular momentum. The above discussions indicate that unstable planets with Oj < 0^2 
expand to overflow their Roche lobe via both the LI and L2 points. The planet's angular 
momentum is absorbed by the lost mass exterior to the orbit and it become tidally disrupted in 
situ. But, unstable planets with a.; > 0,12 overflow their Roche radius via the LI point and they 
must absorb all the initial angular momentum. In this case, the planets' mass loss is terminated 
during stage 3 when they attain at and et- In the limit that et « Cj, the fractional loss of mass 
AMp/Mp ~ {at - Ri)/2aL + ef/{l - ef) 0.05 - 0.3 (see eq. 121). Most short-period planets 
are fractionally less massive than the more distant planets, though this trend is probably due to 
observational selection effects. Nevertheless, the absence of very massive planets among them may 
be indicative of some mass loss for the short-period planets. 

In the above discussions, we show that all planets with P > PeL{e = 1) are stable against 
runaway tidal inflation. Also no planet can continue to overflow its Roche lobe with P > Pile = 1). 
Nevertheless, PL(at,et) can be considerably larger than the circularization period for isolated 
planets with similar kinematic properties. Equation (18) indicates that the decay time of e would 
be extremely short owing to large radii (= Roche radii) during overflow phases: in the case that 
the Roche radius Rl = 4.92i?j at a = 0.04 AU at the beginning of overflow, equation (18) gives Te 
as short as ~ 10""^ years. When the planets' mass loss through Roche lobe overflow is terminated 
at aL, Tg{aL, Rl) is still less than 1 Myr (sec cq. 110). Unstable planets with a > may retain 
their residual eccentricity. Thus, around stars of similar ages, we expect the coexistence of planets 
with periods up to 2-3 weeks, some with initially highly eccentric orbits but which migrated 
outward and therefore currently have circular orbits, and others with orbits in which the planets 
did not inflate too much and have stayed near their original orbits, with initially and currently 
modest eccentricity. 

In the absence of any tidal heating, protoplanets with Rp 2± 2Rj overflow their Roche radius 
when their a is reduced to ~ 0.014 AU. Trifling et al. (1998) suggested that after their tidal 
interaction with their protoplanctary disks brings them to such a distance from their host stars, 
protoplanets lose mass primarily through the LI point. They also proposed that the angular 
momentum transferred from the lost mass to the remaining protoplanet may be adequate to 
compensate for the angular momentum loss to the disk. Since this halting mechanism can only 
occur when the protoplanets overflow their Roche lobe, it is unlikely to be effective for planets with 
periods greater than 3 days, even taking into consideration their inflation by tidal dissipation. For 
protoplanets which migrated to such close range, this mechanism alone cannot sustain a protracted 
angular momentum drainage by a long-lived disk. In the absence of a strong tidal interaction with 
rapidly rotating host stars or the existence of a cavity in the inner regions of their disks to halt 
their migration (Lin et al. 1996), all protoplanets that migrate to the close neighborhood of their 
host stars would perish through tidal disruption. 
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ll! view of the uncertainties in the magnitude of Qp, we briefly discuss the dependence 
of various quantities on it. The quantities which hnearly vary with are the time scale of 
eccentricity damping and tidal inflation. A relatively large Qp can significantly prolong the process 
of tidal inflation. But, the dependence of the critical eccentricity (ex,) for runaway inflation on Qp 
is more modest. For sufficiently large e's, the critical values of a's for the onset of tidal inflation, 
Roche-lobe overflow, and deflation are relatively insensitive to the magnitude of Qp. Since the 
mass reduction factor is determined by the stalling semimajor axis, it too is an insensitive function 
of Qp. Thus, provided the magnitude of Qp < 10^, our overall conclusion that runaway inflation 
may lead to Roche lobe overflow and the disruption or the displacement of ultra short-period 
eccentric extrasolar planets does not sensitively depend on the equilibrium tidal prescription we 
have adopted. 



6. Numerical Simulation of Tidal Inflation 

The analytic approximations presented in the previous sections are useful for categorizing 
various potential evolutionary paths. However, the stability criteria in equations (46) and (49) 
neglect the planets' e damping during their evolution towards a state of thermal equilibrium. The 
discussion in the previous sections also indicates that depending on the planets' initial conditions, 
there are many possible evolutionary paths and that therefore some of the analytic approximations 
may not be appropriate. 

In this section, we present numerical solutions to illustrate the tidal evolution of ultra-short 
period planets. For this quantitative analysis, we consider a simple extrasolar system which 
consists of a solar-mass star and a Jovian planet with an initially eccentric orbit. The equations 
governing the time evolution of e, Rp, and a involve total energy and angular momentum 
conservation of the whole system consisting of a giant planet and a sun-likc host star. Some 
simplification has been made in our approach. We do not perform the calculation for the 
expansion/contraction of a planet that is placed in the Roche lobe potential. Instead, we use 
the evolution code of BLM for a spherically symmetric, gravitationally isolated giant planet to 
estimate the evolution of the planet's radius. We include the heating from both tidal dissipation 
(appropriate for planets with a ranging from 0.02 AU to 0.05 AU) and the stellar irradiation. 
Several evolutionary calculations arc carried out for each value of Mp, for different (constant with 
time) values of dissipation rates. In each case the evolution is calculated up to a time of several 
Gyr or until the radius becomes constant at the value corresponding to thermal equilibrium, that 
is, when the luminosity generated by dissipation equals the luminosity radiated from the interior. 
From these calculations we generate a table that supplies the relation Rp = /(.Etide, Rp, Mp). Then 
we can interpolate/ extrapolate the values of Rp from the table by entering it with any dissipation 
rate -Etide calculated from either the equation of e-damping (eq. 27 with the synchronization 
condition: n = Op), or that of synchronization damping (eq. 27 with the condition for a circular 
orbit: e = 0). We do not run the simulation with both synchronization and circularization 
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heating at the same time. The heating due to synchronization is much smaller than that due 
to circularization, as roughly estimated in equation (44). This point will be justified by our 
simulation results in §6.2. In summary, we solve equations (11), (27), (103), (98), and (126) 
together with the interpolation/extrapolation function Rp = f{Etide,Rp,Mp) for time evolution of 

e, Op, a, Rp, Mp, and Eude- 

Several cautions need to be addressed with our approach. First of all, since the simulation 
carried out by BLM does not take into account the planet spin, we do not include the energy 
Espin to calculate Rp. In §3.1, we already showed that Egpin is smaller than Ei^ind, and therefore 
ignoring Egpin should not influence the internal structure too much. The second caution is that 
the value of Rp in the tables are computed using constant heating rates. In principle, the evolution 
of stratification of planets due to thermal imbalance in the case of quasi-hydrostatic equilibrium 
should react differently to different heating methods. In §3.2, we suggest this assumption might 
be appropriate because much of the gas in giant planet interiors has a highly non-ideal equation 
of state so that the tidal dissipation mainly leads to temperature increases without any significant 
changes in the pressure. 

Mass loss starts to occur once the photosphere of an inflated planet reaches its Roche 
radius. Overflow gas would redistribute mass and angular momentum of the whole star-planet 
binary system and then alter the semi-major axis with time. In our numerical calculations, we 
assume the planets overflow their Roche lobe via the LI point only and that the system evolves 
in a conservative manner. We employ equation (98) to estimate the mass loss rate Mli- Our 
simulations are terminated once Mli > 0. This situation usually occurs when the planet starts 
to contract gravitationally due to inefficient tidal dissipation as e drops to low values in our 
simulations. 



6.1. Eccentricity Damping 

Table 1 shows the minimal e's required for a giant planet to reach its Roche lobe. This e 
is equivalent to in equation (48). While the larger initial planet radii {Rp^init = l-8Rj for 
IMj and Rp^mit = 2.0i?j for 0.63Mj) represent earlier phases of giant planets, the smaller initial 
radius {Rp^init = ^■'^Rj) is used as a later phase of a giant planet that already contracted and is 
reaching thermal equilibrium. In all these models, we set = IMq and Q'p = 10^. The data 
in parentheses are the times taken for the planet to expand from the initial conditions to fill its 
Roche lobe in units of Myr. This time scale is equivalent to tr in equation (45). The table shows 
that this time scale increases with a, resulting partly from larger Tg as indicated in equation (18) 
(i.e. e decays more slowly at larger a due to a weaker tidal effect) and partly from larger Rl at 
larger a (i.e. a planet requires a longer time to reach the larger Roche lobe at larger a). The e's 
for small-size planets are higher than those for large-size planets as a result of larger gravitational 
binding energy for the planets starting from smaller radii. The same reason explains why planets 
with a core require larger e's to reach their Roche radius than planets without a core. The results 
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of our simulations indicate that a planet of 0.63Mj without a core (with a core) at 0.04 AU 
starting at the critical eccentricity 0.140 (0.156) and 2 Rj will reach its Roche radius in about 1 
Myr. This time scale is shorter than the one given by equation (45), which can be explained by 
the steep change of e and Rp as Rp approaches Rl and the reduced gravitational binding energy 
that needs to be overcome for larger radii. 

The results in Table 1 also show the value of the minimum e is inversely proportional to Rp 
which is in agreement with the expression in equation (48). The weak dependence with Mp is 
probably associated with Co being an increasing function of Mp. Although planets located at 
larger a's require larger values of minimal e's and thereby survive more easily, the minimum value 
of e increases less rapidly with a than that expressed in equation (48). The minimum values of e's 
are generally larger than e^,, especially for small a. For a < a^^^ the more appropriate condition 
for planetary inflation \s e > en (see eq. 51). The numerical results are in good agreement with 
the expression of in equation (46). 

Figure 1 illustrates the time evolutions of e, Rp, a, and £'te for a coreless planet with an 
original Mp = O.QSMj, initially located at 0.04 AU (solid lines) and 0.06 AU (dashed lines). 
Starting with e !^ 0.14 (e !^ 0.26), the planet at a = 0.04 AU (a = 0.06 AU) spends t lA Myrs 
{t 8.2 Myrs) and then is able to expand up to its Rl at which = — -E-te- After this point, 
the planet will contract since the dissipation rate Etc has already decreased with time before the 
planet's size reaches its Rl as is shown in the lower-right panel of Figure 1. In the case for larger 
values of a such as 0.06 AU, E'te starts to increase more slowly than that for smaller a as a result 
of a longer damping time scale of eccentricities, leading to a slower rate of expansion as shown in 
the upper-right panel in Figure 1. As long as the initial is large enough for the planet to expand 
slowly through the initial phase and then is able to get larger, the inflation rate is subsequently 
enhanced by its larger radius (£Jte Rp) and smaller gravitational binding energy. The lower-left 
panel in Figure 1 shows that a decreases with time due to circularization. The critical eccentricity 
0.14, calculated by the detailed simulation for a planet of Mp = 0.63Mj at a = 0.04 AU as shown 
in Table 1 and in Figure 1, just coincides with the value br = 0.14 based on the rough estimate as 
shown in equation (46), and is larger than the value cl = 0.06 that is estimated based on equation 
(48) with an extrapolation of BLM's data as described in §3.2. 

Based on the physical process that planet mass is lost through LI (see eq. 98), we calculate 
the smallest initial e^'s required for a giant planet to move from its initial orbit with a semi-major 
axis ainit to its final orbit with a semi- major axis a final- The results without a core and with a 
core are shown in Tables 2, 3, and 4. The mass loss rate is calculated based on equation(98). 
The data for A = 1 and A = 3 are very close to the data for A = 2. Therefore we only show the 
results for A = 2 in this paper. Giant planets move away from their host stars by receiving orbital 
angular momentum transfered to them from the overflow gas, and their masses are reduced by a 
fraction before they become detached from their Roche lobes. The data in bold face as displayed 
in parentheses show the final masses of giant planets in units of IMj after their radii become 
detached from their Roche radius. Typically, Mp is reduced by « 15% to 40% for planets with 
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initial e's comparable to the critical values. The surviving mass is expected to be smaller if the 
planets' initial e's are larger than the critical values. These minimal e's depend on the initial 
conditions: the greater the initial gravitational binding energy (due to smaller radii, larger mass, 
or the existence of a core), the greater the dissipation, and therefore the larger the e, is needed to 
inflate the planet to attain the same final condition. For an outward migration from 0.03 AU to 
0.04 AU, a planet starting with IMj and 1.2Rj requires e ^ 0.2, a planet initially with IMj and 
l.SRj needs e ^ 0.17, and a planet with 0.63Mj and 2Rj should start with e ^ 0.14. However, a 
giant planet of IMj and l.SRj originally located at 0.04 AU requires e > 0.204 (without a core) 
or e > 0.223 (with a core) to migrate to 0.05 AU, indicating the trend of difficulty for a giant 
planet located further away to move outwards by Roche lobe overflow. The difference of critical e's 
between a planet with a core and without a core increases with a because weaker tidal interaction 
occurs at larger a. The first values in parentheses (not in bold face) show the time scale in units 
of Myrs for a from ttinu to a final- The time scales range from some small fraction of one Myr to 
more than 10 Myrs depending on initial values of e, ainu, and a final- 

Figure 2 displays the time evolution of e, Rp, a, Mp, and Ete for a planet without a solid 
core. The initial radius is 1.8 Rj and the initial mass is 1 Mj. Starting with the e = 0.155 which 
is the minimal eccentricity to get to the final orbital distance 0.04 AU as shown in Table 2, the 
planet reaches its Roche lobe around t = 0.34 Myrs, marked by dashed lines in the figure. The 
beginning of the Roche overflow phase is illustrated in the Mp/Mj vs t plot where Mp starts to 
fall off around t = 0.34 Myrs. At the same time, Ete reaches the maximal value. Subsequently, 
the planet loses its mass, moves outwards, and increases its size as a result of Roche lobe overflow. 
The decrease in Ete during the overflow phase is a consequence of the increase in a. The mass loss 
rate is determined by equation (98) with the condition that A = 2. After having lost 0.16 Mj, 
the planet starts to contract and therefore detaches from its Roche lobe at a « 0.04 AU due to 
inefficient Ete (and e) when t pa 0.45 Myrs. Note that a slightly decreases with time before the 
beginning of the Roche phase {t ^ 0.34 Myrs) as a result of orbital circularization. However, the 
change of a is dominated by the effect of overflow during the mass- losing phase (i.e. Id^l > |de| in 
equation (96) when t > 0.34 Myrs). 

Figure 3 illustrates the evolution of the same planet as shown in Figure 2 but with a larger 
value of initial eccentricity. With a large value of initial e, the whole dynamical evolution should 
more or less go through a sequence of the stages as illustrated in §5. Owing to the larger initial 
e = 0.2 shown in Figure 3, the planet loses more mass and therefore migrates even farther out than 
the case (as shown in Figure 2) with the initial e = 0.155. The planet enters its Roche-overflowing 
phase at t pa 0.13 Myrs. During the period t = 0.13 to 0.2 Myrs {i.e. the early phase of the Roche 
overflow), Rp increases roughly from 4Rj to more than 7Rj but e declines approximately from 
0.14 to 0.12, demonstrating "Stage 2" described in §5. The following phase after t > 0.2 Myrs 
featured by the flattening trend of Rp with time confirms the stagnation stage (i.e. Stage "3") 
described in §5. The final contraction stage is not able to be simulated with our method which is 
based on extrapolation/interpolation of the planet data inflated up to the equilibrium radii. 
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The above results are based on the assumption that overflow occurs only through the LI 
point. However, there might be an additional concern in the star-planct interacting system. The 
gravitational potential difference between the LI and L2 points is drastically reduced for the 
small mass ratio Mp/M* ^ lO'^, therefore raising the possibility of overflow through the outer 
Lagrangian point (L2). Our simulation for a planet of IMj without a core at a = 0.03 AU shows 
that for the initial eccentricity e = 0.155 (e = 0.5), the rate of mass transfer across equipotential 
surfaces is about Ma ~ 4 x 10^^ g/s (9 x 10^^ g/s), which should be roughly comparable to the 
mass flux via the LI point. This flux is given by equation (81). Let us estimate this mass flux Mn 
based on the conditions at the photosphere: p = 10^^ g/cm'^, Cg = 4 x 10^ cm/s. These values 
would give rise to the mass loss rate Mli ~ 10^^ g/s. Hence M^i ^ Ma for initial e = 0.155, but 
Ma ^ Mli for initial e = 0.5. As we discuss in §4, mass should be primarily transferred via the LI 
point in the small-e limit since the bottom of the planet's photosphere underfills the LI point and 
the density scale height hp is smaller than AD as shown in equation (66). In the large-e regime, 
the bottom of the planet's photosphere overfills LI, probably leading to a significant amount mass 
loss through L2. If the planet starts at a = 0.02 AU, our simulation shows that Ma ^ 5 x 10^^ 
g/cm^ for initial e = 0.25, therefore increasing the chance of overfilling the LI point as indicated 
by equation(84) for smaller a. 

In the previous section, we suggested that in a system with 2ea < hp, mass transfer is 
continuous because the gas removed from the atmosphere is replenished on a dynamical time scale 
by the requirement of hydrostatic equilibrium. For planets with larger e's, the envelope as well as 
the atmosphere needs to adjust. Also as reasoned in the previous section, the underlying envelope 
of the planet expands at a rather slow rate, which should not be of importance in determining M 
if there is continuous mass overflow in one orbital period. Our simulation shows that the fastest 
expansion rate of the planets' envelope is Rp ~ 0.1 cm/s for an initial e = 0.5, and ~ 0.01 cm/s 
for an initial e = 0.155 in the case of a giant planet of IMj without a core at a = 0.03 AU. 

When a planet moves from perihelion to aphelion, its Roche radius increases at a rate 

which is indeed much larger than the expansion rate ~ 0.1 cm/s due to the e-damping from our 
simulation. On the other hand, 

2ae Ri 1.66 x 10^ ( — ) ( ) cm - /i. 2 - 3 x 10^ cm, (124) 

VO.2; V0.04AU; ' ^ ^ 

where hp is evaluated at 0.04 AU (see §4.2). As a result, one would expect mass overflow via 
the LI point to be sustained in different phases of one orbital period despite the re-expansion of 
planet's atmosphere, although the most prominent overflow occurs around the perihelion. Hence 
we conclude that using the Roche radius based on perihelion is a reasonable approximation for 
small e's. However, hp would be much larger than 2ae at smaller a since hp gets larger due to 
more intensive stellar radiation and furthermore 2ae gets smaller. 
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So far we have not considered the effect of the planet's tidal interaction with the protostellar 
disk and the consequent inward migration. The primary concern about inward migration is that if 
the planets' orbits are almost circularized by the time they reach a = 0.04 AU, tidal dissipation 
within them wouldn't be available to inflate their radii. As has been discussed in §3.3, there might 
be a couple of theoretical mechanisms to excite/damp e as giant planets migrate in through a 
disk. We do not include these uncertain mechanisms in this paper, but only take into account the 
damping of e due to as a function of various inward migration rates, using a simple prescription 
a/ a =constant. Figure 4 shows the evolution of e of a corelcss synchronized planet, with IMj 
and 1.8i?j, moving toward its host star from a = 0.4 AU on three different time scales: —a/ a =5, 
1, and 0.1 Myr. The simulations start with e = 0.2 and e = 0.3 and end before the planet reaches 
its Roche lobe. As shown in Figure 4, e's do not drop noticeably from their initial values until 
the planet reaches a = 0.04 AU; thus at this point enough energy from orbital circularization still 
remains to inflate the planet to its Roche radius. The orbital distance a for Roche overflow shifts 
to smaller values when a/a oi initial e is smaller. 

The minimal e's listed in the first column (i.e. those data with the initial condition ainit = 0.02 
AU) in Table 4 are unexpectedly large compared with those in the first columns in Tables 2 and 3 
when the values of these minimal e's are ^ 0.2. In principle, an internally-heated planet starting 
with a smaller mass and a larger size should overcome its gravitational binding energy more easily. 
In reality, this surprising outcome results from a special initial condition for any large e: the 
initial radius Rp = 2Rj is already larger than the planet's Roche radius Rl for Rp = 0.63Mj, 
a = 0.02 AU, and e > 0.2 according to equation (34). Under such a circumstance, the planet 
cannot be inflated to a large size before overflow happens, leading to less intensive tidal heating 
and therefore requiring a larger e to reach a = a final- In contrast, a planet at ainit = 0.02 AU 
starting with more mass such as IMj can be inflated to a larger size before overflow starts to take 
place, resulting in more intensive tidal heating and thereby requiring a smaller e to migrate to the 
same final location a final- 

6.2. Synchronization 

In this subsection, we shall show that unlike orbital circularization, synchronization is not 

a promising source of dissipation for a planet within 0.04 AU of its host star since the planet 
is almost completely synchronized during the course of its inward migration. In this analysis, 
wc consider the dissipation solely due to synchronization, but including orbital migration and 
the planet's inflation/contraction. There will be three equations governing the time evolution 
of Qp, Rp, and a. We employ the same assumptions and approach for evaluating Rp and a for 
synchronization damping as we described in the previous section for eccentricity damping. In the 
case of synchronization for a circular orbit, the dissipation is given by E\,q in equation (28): 
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Lq. (125) 

In equation (5), Op contains two pieces, each of which is associated with a different physical 
process: synchronization based on the conservation of orbital and spin angular momentum (e.g. 
see Murray & Dermott 1999), and inflation/contraction based on the conservation of spin angular 
momentum such that 

% = ^tide + ^|/pO=constant = ( sign(n - ^p)— - Yl ^P' '^"^^^^ 

\ TCI IpJ 

Figure 5 shows the time evolution of a, Op, rp, ErU, and the synchronization factor 0,p/n, for 
a 0.63 Mj planet with three different migration rates a/a = 5, 1, and 0.1 Myr. Initially the planets 
are fast rotators at a spin period 27r/Op = 2 days, and they start from 0.4 AU with Rp = 2Rj. 
This corresponds to the condition Op/n = 46.48. Then they migrate until they overflow their 
Roche radii. 

In the beginning the spin periods 27r/Op decrease as the planets move in as a result of 
contraction. The contraction phase continues until the tidal interaction becomes so important 
that the spin frequencies Op are slowed down by the orbital frequencies n which increase with 
time due to migration. After n catches up with Op (this occurs around 0.05-0.1 AU depending 
on migration rates) and then becomes larger than Op as a result of continuous migration, the 
spin periods decrease again mainly due to the fact that Op tries to catch up with n because of 
synchronization. As a result there are two dissipation peaks due to synchronization during the 
course of migration as illustrated on the Eiq/Lq diagram in Figure 5. The first dissipation peak 
happens at a relatively large o, when n catches up with Op, and the second peak occurs at a 
smaller a when Op catches up with n. The second peak has a relatively small amplitude compared 
to the first peak because the difference between n and Op is small during the second dissipation 
peak. Therefore, one can essentially state that the system has reached its synchronous state after 
n catches up with Op at around 0.05 AU. This result is in agreement with the discussion in §3.2. 

The continuous contraction as shown in the upper-right panel in Figure 5 is just a result of 
small Ttn shown in the lower-left panel in Figure 5. A fast migration rate a/a =0.1 Myr is unable 
to boost tidal dissipation high enough to inflate the planet. The final Roche lobe overflow taking 
place at a ^0.015 AU (0.016 AU for -a/a = 0.1 Myr, 0.015 AU for -a/a =1 Myr, and 0.014 
AU for —a/ a = 5 Myr) is therefore not due to radius expansion, but due to a decrease of the 
Roche radius as the planet migrates in. Although the simulation is carried out only for a coreless 
planet with a relatively small mass Mp = 0.63Mj, planets with cores or of larger masses such as 
Mp = IMj are more difficult to inflate to their Roche radii by synchronous dissipation due to 
their larger gravitational binding energy. 

We also plot the evolution of a coreless planet with 0.63Mj which is being tidally heated as 
a result of synchronization, starting without any orbital migration. The result is illustrated in 
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Figure 6. The planet starts as a rapid rotator with Q,p/n pa 7.64 at a = 0.03 AU (i.e. 2-!t/^p = 0.25 
days). This value of i}p has been chosen to be sufficiently large for the generation of a substantial 
amount of EfQ, but small enough to limit -Egpin to be less than -Ebind- The synchronization time 
scale is extremely short: tq ~ 1000 yrs according to equation (17) for Mp = 0.63Mj, a = 0.03 
AU, Rp = 2Rj, and 27r/Qp = 0.25 days. Figure 6 shows that Up/n — > 1 within a few thousand 
years. Consequently, even though the equilibrium radius for the intense initial dissipation rate 
(O.OIL0) is Rp = 2.QRj, the planet exhausts its energy source and stops its expansion in a couple 
of thousand years, well before reaching its Roche lobe. This result confirms the general statement 
TRn > Tfi as described by equation (44). 



7. Summary and Discussion 

Based upon interpolation/extrapolation of numerical results of internal structure by BLM, 
we have demonstrated, analytically and numerically, that giant planets with masses < IMj, 
initial radius l.SRj, modest e's (;^ 0.1 — 0.22), and small a's (< 0.04 AU) are likely to be 
tidally inflated beyond their Roche lobe as they undergo e damping. These conditions may be 
slightly modified by their internal structure (with or without cores) and by inward migration due 
to disk-planet interaction. The subsequent mass loss depends on the modification of the internal 
structure after these tidally infiatcd giant planets have lost some of their initial mass. Therefore 
this adjustment is difficult to estimate through our approach based on the numerical data without 
mass loss by BLM. The degree of outward migration will determine whether or not this Roche 
lobe overflow model is likely to be responsible for the lack of giant planets within 0.04 AU. How 
much mass the planet will lose and how far out the planet will migrate is a matter of competition 
between several factors. Mass loss increases the expansion rate for a given dissipation rate. Our 
simulation shows that the planet radius, in this case confined by the Roche radius, increases as the 
planet gains specific angular momentum from the overflow gas through the LI point. Increasing 
the semi-major axis a due to Roche overflow via the LI point inevitably reduces the internal 
tidal heating. All of these factors undoubtedly depend on the evolution of the internal structure 
under the condition of mass loss. Nevertheless, we adopt a simplified method by employing the 
parameter A to relate the mass loss rate to the expansion rate of the planet (see cq.(98)). Based 
on reasonable values of A between 1 and 3, our results indicate that a planet located at larger ainit 
requires a larger value of critical e to move out to the same a final- The minimal values of e to 
reach the Roche radius increase rapidly as ainit increases from 0.03 AU to 0.05 AU as a natural 
result of the steep dependence on a for tidal interaction. This result depends very weakly on 
the magnitude of the highly uncertain value of Q'p. Therefore, our model suggests that the lack 
of planets inside 0.04 AU might be explained by outward migration due to Roche lobe overfiow 
rather than extreme mass loss. On the other hand, if the protoplanetary disk is still present, the 
outward migration could be countered by the inward migration from disk-planet interaction; in 
this case extreme mass loss could actually provide the explanation. The eccentric orbits e > 0.2 of 
extrasolar Jovian planets are theoretically likely as described in §3.3 and are commonly observed. 



-44- 



except near the host star. Also, we have shown that eccentricity damping before a planet migrates 
to within 0.04 AU is so small that the mechanism of eccentricity damping can provide sufficient 
energy to inflate a planet to the Roche radius inside 0.04 AU of its host star. Thus Roche lobe 
overflow due to eccentricity damping seems to be a reasonable mechanism to explain the apparent 
lack of giant planets within 0.04 AU based on our approach. Around stars of similar ages, we 
expect the coexistence of planets with periods up to 2-3 weeks, some with initially highly eccentric 
orbits but which migrated outward and therefore currently have circular orbits, and others with 
orbits in which the planets did not inflate too much and have stayed near their original orbits, 
with initially and currently modest eccentricity. 

The amount of mass loss provides an estimate of the metallicity enhancement of the host 
star, provided that (1) the planet's envelope is substantially enhanced in metals relative to the 
stellar metal abundance, and (2) some fraction of the overflow mass is able to accrete onto the 
star. The internal structure of Jupiter and Saturn is somewhat uncertain (Wuchterl et al. 2000) 
ranging from that with a solid core which contains most of the heavy elements to those where 
metal-enriched gas is well mixed in the convective envelope. But in both sets of models, the heavy 
element abundance inside Jupiter and Saturn is thought to be 5-10 times larger than that of the 
Sun. During the stage 2, if the tidally inflated planets lose most of their mass via the LI point, all 
the tidal debris would be accreted by the host stars. But if mass flow occurs through both the LI 
and L2 points, the tidally inflated planet would be disrupted rapidly with the tidal debris forming 
an accretion disk which eventually channels the metal rich gas onto the host star. As long as the 
stars do not have a deep convective envelope, as is true in the case of early-type stars or during 
the later evolutionary stages of solar-type stars, the effect of metallicity enhancement should be 
observable on the stellar surface (Sandquist et al. 2002). Further investigation should be carried 
out to clarify these effects. 

Once giant planets inwardly migrate to ~ 0.04 AU through the planet-disk interaction, their 
spin frequencies are more or less synchronized with orbital frequencies. In spite of the fact that 
some dissipation is generated as spin frequencies catch up with orbital frequencies within 0.04 AU 
as a result of inward migration, our results with different inward migration rates show that the 
internal tidal damping due to synchronization alone is weak enough that the giant planets contract 
during their entire evolution. A simulation for a fast-spinning planet at 0.03 AU {0,p/n ~ 7.64) 
with no inward migration due to disk-planet interaction shows that the giant planet still cannot be 
inflated to its Roche radius, consistent with the relation that tr^ > tq induced from equation (44). 
However, planet winds may enhance synchronization heating if they can carry a sufficient amount 
of spin angular momentum out of the planet, thereby holding the planet out of synchronous 
rotation. This effect should be considered in the future, in a model perhaps analogous to that of 
the black widow pulsar (Applegate & Shaham 1994). 

The migration model through the planet-disk interaction alone seems to have difficulty in 

accounting for an over-abundance of lower mass planets at a < 0.1 AU (Armitage et al. 2002). 
Although the over-abundance may possibly be a result of observational bias, it could, particularly 
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near a = 0.04 AU, be due to the retreat of Roche-lobe-filling planets from the region inside 
a = 0.04 AU. Whether the paucity of the massive close-in Jovian planets (Zucker &: Mazeh 2002) 
results from our scenario is left undetermined until more simulations are carried out for the 
close-in planets with Mp > IMj (but, see Paetzold & Rauer 2002 and Jiang et al. 2003for another 
model). The distribution of giant planets at 0.03 AU has been simulated by Trilling et al. (1998) 
by taking into account the effect of planet-disk interaction, the tidal torque from a fast rotating 
star, and Roche lobe overflow, but without considering tidal inflation. The effect of tidal inflation 
causes giant planets to reach the Roche lobe phase at a larger a. 

The expansion/contraction rates Rp/Rp used in our simulation are taken from the simulations 
by BLM which are carried out for constant heating rates. But, the dissipation rates actually 
increase with time for e larger than the minimal values as a consequence of expansion rates larger 
than eccentric damping rates: 



We have used equations (29) and (18) to derive the above equation. Therefore, the minimal 
values of eccentricities required for planets to reach the Roche limit in our simulations are likely 
underestimated. However, one thing to keep in mind is that all the uncertainty of tidal friction is 
hidden in the dissipation quantity Qp, and we have assumed that the tidal heating is uniformly 
distributed in planet interiors. Since tidally inflated planets have radiative envelopes (BLM) and 
it is the radiative envelopes that respond most to the deposition of thermal energy, the effect of 
tidal expansion might be enhanced due to more dissipation generated in the radiative envelopes 
as a consequence of the action of radiative damping on the dynamical tide (Zahn 1977). 

The overflow mechanism that we explore in this paper might offer a reasonable explanation 
for the recent HST observation results on a lack of Jovian planets with orbital periods ^ 4 
days around main-sequence stars in the globular cluster 47 Tucanac (Gilliland ct al. 2001). In 
addition to the breakup of planetary systems during frequent encounters with single stars and 
binaries in the dense core of 47 Tuc or the scenario that the lack of giant planets results from 
metal-poor circumstellar disks, eccentricities of planets can get excited after a close encounter 
(Davies & Sigurdsson 2001) with a binary system. Our simulations show that the minimal 
eccentricities required for a giant planet of IMj starting with \.2Rj to move from ainit = 0.04 
AU to a final = 0.05 AU is ~ 0.285, and that the minimal eccentricity for the same planet but at 
o-init = 0.05 AU to reach its the Roche radius is ~ 0.3. Therefore, with large eccentricities possibly 
excited in dense stellar environments, a lack of the Jovian planets with very short periods in 47 
Tuc could result from the same effect via Roche lobe overflow caused by sufficient tidal infiation. 
Further investigation on this topic is ongoing (Oh & Lin 2002). Tidally- infiated Jovian planets 
in eccentric orbits should produce deep transit light curves with shorter flat bottoms, providing 
an unique signature for the surveys for transiting hot Jovian planets in young extrasolar systems 
with spectroscopic follow-ups. However, the possibility of finding these objects in surveys, even 
conducted in some young open clusters, would be small due to the short time scales of tidal 



dEte/dt 
Etc 




(127) 
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infiation in our model. 
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Table 1. Minimal eccentricities required to reach Roche lobes. The cases without a core (the 
values before /) and with a core (the values after /) are displayed. The numbers in parentheses 
are the time taken from the initial conditions to the initial Roche phase. 



Mp 


Rp,init 


0.02 AU 


0.03 AU 


0.04 AU 


0.05 AU 


IMj 
IMj 
0.63Mj 


l.SRj 
1.2Rj 
2.0Rj 


0.095(0.08)/0.097(0.07) 
0.129(0.20)/0.129(0.18) 
0.055(0.05)/0.055(0.06) 


0.147(0.49)/0.154(0.41) 
0.183(1.64)/0.192(1.07) 
0.101(0.32)/0.103(0.23) 


0.203(1.89)/0.220(1.61) 
0.252(7.17)/0.284(7.18) 
0.140(1.35)/0.156(0.93) 


0.274(4.44)/0.299(3.51) 
0.313(12.7)/0.385(12.6) 
0.195(3.82)/0.228(2.65) 



Table 2. Minimal eccentricities for planets of initial IMj and l.SRj to migrate outwards from 
initial separations ainit to final separations a final based on the equation (93) with the condition 
that A = 2. The first values in parentheses show the time scales in units of Myrs for a to reach 

agnai from flinit- The second values in parentheses, which arc in bold face, show the final masses in 
units of IMj when planet radii detach from their Roche lobe. The cases without a core (the 
values before /) and with a core (the values after /) are displayed. 





Oinit = 0.02 AU 


ttinit = 0.03 AU 


ttinit = 0.04 AU 


>0.04 AU 

>0.05 AU 
>0.06 AU 


0.131(0.19, 0.70)/0.147(0.15, 0.70) 
0.145(0.29, 0.62)/0.170(0.20, 0.62) 
0.152(0.59, 0.57)/0.183(0.47, 0.57) 


0.155(0.45, 0.84)/0.167(0.36, 0.85) 
0.162(0.45, 0.76)/0.180(0.34, 0.76) 
0.174(0.47, 0.69)/0.199(0.32, 0.69) 


0.204(1.83, 0.87)/0.223(1.45, 0.87) 
0.207(1.69, 0.80)/0.229(1.24, 0.79) 


Table 3. 


Same as Table 2 for a planet of IMj, but with a different initial radius 1.2Rj. The 


eccentricities shown here are larger than those in Table 2 owing to the fact that a planet of a 




smaller size has larg 


;er gravitational binding energy. 




flfinal 


ttinit = 0.02 AU 


ttinit = 0.03 AU 


ttinit = 0.04 AU 


>0.04 AU 
>0.05 AU 
>0.06 AU 


0.157(0.26, 0.69)/0.167(0.20, 0.70) 
0.168(0.36, 0.62)/0.188(0.28, 0.62) 
0.174(0.65, 0.56)/0.199(0.50, 0.57) 


0.189(1.44, O.84)/0.200(1.31, 0.85) 
0.195(1.31, 0.75)/0.209(1.05, 0.76) 
0.203(1.18, 0.69)/0.223(0.80, 0.69) 


0.253(7.00, 0.84)/0.285(6.62, 0.81) 
0.255(6.68, 0.79)/0.286(6.17, 0.78) 
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Table 4. Same as Tabic 2, but with a different initial planet mass 0.63Mj and a different initial 
radius 2Rj. The eccentricities shown here are in general smaller than those listed in Table 2 and 
Table 3 except for the cases starting with ainit = 0.02 AU. 



flfinal 


Oinit = 0.02 AU 


ttinit = 0.03 AU 


ainit = 0.04 AU 


>0.04 AU 

>0.05 AU 
>0.06 AU 


0.124(0.58, 0.44)/0.136(0.27, 0.44) 

0.145(0.87, 0.39)/0.177(0.33, 0.39) 
0.171(1.16, 0.36)/0.232(0.36, 0.36) 


0.122(0.42, 0.54)/0.137(0.29, 0.54) 

0.135(0.50, 0.49)/0.160(0.33, 0.48) 
0.153(0.31, 0.44)/0.194(0.36, 0.44) 


0.150(1.18, 0.55)/0.173(0.71, 0.56) 
0.159(1.08, 0.51)/0.190(0.60, 0.51) 
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Fig. 1. — Time evolutions of e, Rp, a, and Ete for a planet of a original Mp = 0.63Mj without 
a solid core initially located at 0.04 AU (solid lines) and 0.06 AU (dashed lines). Starting with 
e « 0.14 (e =f« 0.26), the planet at a = 0.04 AU (a = 0.06 AU) spends t 1.4 Myrs {t 8.2 Myrs) 
and then is able to expand up to its equilibrium radius which is just equal to its Rl- 
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Fig. 2. — Time evolution of e, Rp, a, Mp, and Etc for a planet with original Mp = IMj without a 
solid core initially located at 0.03 AU. The mass loss is calculated by using equation (98) with the 
condition that A = 2. The moment t = 0.34 Myrs marked by a vertical dashed line indicates the 
beginning of the Roche overflow phase. 
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t (Myrs) t (Myrs) 

Fig. 3. — The same simulation as the one shown in Fig. 2 except starting with a higher value of 
e = 0.2. The moment t = 0.13 Myrs marked by a vertical dashed line indicates the beginning of 
the Roche overflow phase in this case. 
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Fig. 4. — Evolution of a for a planet of 1 Mj and 1.8 Rj without a core and starting at a = 0.4 
AU, moving inwards at three migration speeds: —a/ a = 5 Myr (dotted lines), 1 Myr (dashed lines), 
and 0.1 Myr (solid lines). The simulations start with e = 0.2 and e = 0.3 and end before planets 
reach their Roche lobes. 
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Fig. 5. — Time evolution of Op, Rp, Em, and Q.p/n for a 0.63 Mj planet with no core starting from 
0.4 AU with three different migration rates: —a/ a =5 Myr (dotted lines), 1 Myr (dashed lines), 0.1 
Myr (solid lines). 
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Fig. 6. — Time evolution of $7p, Rp, Etn, and Qp/n for a 0.63 nij planet with no core located at 
0.03 AU with no migration. The planet starts out as a fast rotator with the ratio Qp/n ~ 7.64, and 
then is rapidly synchronized {^p/n — > 1) by the tidal effect within several thousand years. 



